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ABSTRACT 

We have observed the nucleus of the nearby luminous infrared galaxy NGC 4418 with subarcsec 
resolution at 860 and 450 fj,m for the first time to characterize its hidden power source. A ~20 pc 
(O'.'l) hot dusty core was found inside a 100 pc scale concentration of molecular gas at the galactic 
center. The 860 /im continuum core has a deconvolved (peak) brightness temperature of 120-210 K. 
The CO (3-2) peak brightness temperature there is as high as 90 K at 50 pc resolution. The core has 
a bolometric luminosity of about 10^^ Lq, which accounts for most of the galaxy luminosity. It is 
Compton thick (iVn > 10^^ cm~^) and has a high luminosity-to-mass ratio ^500 L© Mq~^ as well as 
a high luminosity surface density lo8-^='=°-^ Lq pc~^. These parameters are consistent with an AGN 
to be the main luminosity source (with an Eddington ratio about 0.3) while they can be also due to 
a young starburst near its maximum L/M. We also found an optical color (reddening) feature that 
we attribute to an outflow cone emanating from the nucleus. The hidden hot nucleus thus shows 
evidence of both an inflow, previously seen with absorption lines, and the new outflow reported here 
in a different direction. The nucleus must be rapidly evolving with these gas flows. 
Subject headings: galaxies: active — galaxies: evolution — galaxies: individual (NGC 4418) — 
galaxies: ISM 



1. INTRODUCTION 

There are galaxies hiding a compact luminous source 
in their nuclei and NGC 4418 is their local prototype. 
The extremely extinguished power source i n this early- 
type disk galaxy was first recognized by iRoche et al.l 
(|1986[) who singled out the galaxy for its infrared bright- 
ness. They found very deep silicate absorption at 9.7 
/im toward the infrared-bright nucleus and deduced that 
the source of the infrared luminosity {^IO^^Lq) is com- 
pletely obscured by a dust shroud of Ay » 50 mag. They 
further inferred from the weakness of line features in the 
infrared and optical spectra that the luminosity source is 
an accreting massive black hole (i.e., active galactic nu- 
cleus or AGN), although they also mentioned a nascent 
starburst embedded within a dusty cocoon as a possible 
alternative. 

Subsequent observations to be reviewed in HI. II have 
supported the presence of a deeply buried compact lumi- 
nosity source, although its nature is still debated. NGC 
4418 is the most absorption-domi nated galaxy in th e 
mid-infrared diagnostic diagram of lSnoon et all ([2001 . 
Its location on the diagonal branch in the diagram sug- 
gests that the luminosity source is almost completely cov- 
ered by opaque dust. Usual spectroscopic signs of AGN 
can be completely extinguished by this heavy absorption. 
Hence, the main piece of evidence for a hidden AGN, 
among other indirect ones, has been the compactness of 
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the source inferred from the absorption depth. It has 
been argued that dust is unlikely to cover an extended 
starburst s o completely as to cause the deep absorption 
featu res (Roche ct al. ' 19861: iDudlev fc Wvnn-Williami 
Il997t |Spoon_ei_ah 2Q0j|TOn the other hand, the ar- 
guments for very young starburst are based on the low 
radio to far infrared ratio and the detection of a warm 
molecular gas absorber toward the nucleus. The former 
ratio is lower than in radio-quiet AGNs and is consis- 
tent with a sta.rburs t in the pre-supernova age (< 5 Myr) 
(jRoussel et al.l l2003^ . The absorbing molecules appear to 
be in the region where X-rays wo uld have destroyed t hem 
if the heating is due to an AGN (jLahuis et al.ll2007[ ). 

The true nature of the luminous nucleus is important 
because this type of heavily exti nguished nuclei could 
double the number of local AGN ()Maiolino et al.ll2003D 
or it could represent a young phase of unusually in- 
tense and compact starburst (^ 100 pc). Luminous 
hidden nuclei of this kind are also important because 
they may be in a short transition phase to unobscured 
AGNs, i.e., in the process of removing the obs curing 
gas ( e .g., through an AGN -driven outflow, Sandcrs 'et al.l 
119881 iFischer et ani2010t ). Such feedback is regarded as 
a key process in galaxy evol ution and hence requires de- 
tailed study (|Hopkins et al.l 120061 : IWeiner et al.1120091 ). 

A new approach to characterize and diagnose the com- 
pact, luminous, and extinguished nuclei of this kind is 
through submillimeter-wave observations at high spatial 
resolution. Sub-mm observations directly probe the ther- 
mal dust emission that constitutes most of the bolomet- 
ric luminosity of such a nucleus. Unlike optical and in- 
frared observations, one can look deep into the obscur- 
ing material since dust opacity decreases at longer wave- 
lengths. Radio interferometry is available in the sub- 
mm to achieve sub-arcsec resolution that is currently un- 
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available in far-IR where the radiation of such a nucleus 
peaks. It is possible to estimate the bolometric lumi- 
nosity and luminosity density of the nucleus using the 
Stefan-Boltzmann law by measuring the source size and 
brightness temperature at a frequency where the nucleus 
has a photosphere. One can also probe the nucleus by 
measuring its mass from rotation of circumnuclear molec- 
ular gas, which can be observed using sub-mm molecular 
lines. The luminosity-to-mass ratio, iboi/A^dyn, so ob- 
tained will further constrain the nature of the nucleus 
since an AGN can have a higher L/M than stars. Molec- 
ular lines can provide additional information on temper- 
ature, density, chemical composition of the interstellar 
medium (ISM), and on the radial motion of the gas af- 
fecting the evolution of the dusty shroud. We had applied 
this approach to the ultraluminous infrared galaxy Arp 
220, which shares some of the above mentioned charac- 
teristics of NGC 4418 including the deep silicate absorp- 
tion. We found that the brighter one of the twin nuclei 
of Arp 220 has Lboi> 2 x 10 " Lp, and Lboi/Mdvn > 400 
Lq I Mq in the central 80 pc (jSakamoto et al.|[2008l here- 
after S08) and that each nucleus has a molecul ar out- 
flow c ausing molecular P-Cygni profiles (Sakamoto et al.l 
[20091) . 

In this paper, we apply the same approach to NGC 
4418 to quantitatively constrain the nature of the deeply 
extinguished nucleus through sub-mm high-resolution 
observations. We made the first sub-arcsecond resolu- 
tion imaging of the galaxy in the submillimeter band, ob- 
serving 860 /im continuum and molecular lines including 
CO (3-2), HCN(4-3), and HCO+(4-3) and also imaging 
the nucleus with 450 //m continuum. We also analyzed 
optical multi-color images to probe circumnuclear dust 
distribution for signs of disturbance and feedback. The 
obtained parameters of the nucleus are similar to the ones 
in Arp 220. They are consistent with any of the follow- 
ing being energetically dominant: an AGN, a young and 
compact starburst, and an AGN-compact starburst com- 
bination. Our results support the previous studies using 
an independent method and indicate a possible path to 
break the degeneracy about the luminosity source with 
future observations of higher quality. We also found a 
sign of outflow from the nucleus in the optical data. In 
the following, we briefly review properties of NGC 4418 
and previous studies on the hidden nucleus in fJlTTJ We 
describe our SMA observations in Sj2] and present the re- 
sults in Sj3l where we also report our examination of op- 
tical images of the galaxy. The results are discussed in 
21 and the conclusions are summarized in [JSJ Our com- 
panion paper reports high-resolu tion 1 mm and cm- wave 
observations of the same galaxy ([Costagliola et aLll2012l 
hereafter C12). 

1.1. NGC 4418 

Stud ies on NGC 4418 and its nucleus since lRoche et al.l 
(jl986[) are summarized below. (To guide readers, an 
optical image and key parameters of the galaxy are in 
Fig. [1] and Table (U respectively. We adopt the galaxy 
distance of 34 Mpc at which 1" is 165 pc.) (1) A 
high infrared-to-optical luminosity ra tio Lir/Lb ~ 20 
(jArmus et al.lll987l: lYoung et aIlll989D and the shape of 
the spectral energy distribution (SED) suggest that most 
of the luminosity of the galaxy is absorbed and rera- 
diated by dust as thermal emission ([Roche fc Chandleil 



[1991 . (2) NGC 4418 has warm mid-IR and far-IR colors 
5'(25/ im)/5'(60/xm) = 0.2 20 and S{60fim)/Sil00fim) = 
1.37 ([Sanders et al.l [20031) . The former puts the galaxy 
in the 'warm' popula tion that often host active nuclei 
(|de Grijp et al.l ll985t iSanders et al[ll988| ). The latter 
also suggests warm dust (color temperature of ^^50 K) 
indicative of reprocessed dust emission. (3) The deep ab- 
sorption has been confirmed by further spectroscopy and 
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100 mag (Roche et al.' '1986'; 'Dudlcv & Wvnn-Wilhami 
1997; Spoon ct al. 2001; Sicbcnmorgcn ct al. 200^. A 
radial gradient of temperature in the nucleus may 
also contri bute to the deep lOum depression in the 
same way iKwan fc Scovillel (I1976D mo deled for proto- 
stars (iDudlev fc Wvnn-Williamslll997[ ). It is also sug- 
gested from a shallow SED slope in the FIR that the 
dust has a high enough column den sity to be optically 
thick even at or beyond 100 /im ([Roche &: C handlen 
Il993t iLisenfeld et"alTl2000[ ). Likely related to this is 
an unusually low ratio of the [C Il j 158 fim line 
to far-IR luminosity in NGC 4418 ([Malhotra et al.l 
Il997t iGracia-Carpio et al.l 120 lit ). (4) The compact- 
ness of the infrared nucleus has been directly con- 
firmed with mid-IR imaging at subarc s econd reso- 
lutioiis ([Wvnn-WiUiams fc BecklinI 119931: lEvans et ahl 
l2003t iSiebenmorgen et a l. 2008). According to the Keck 
imaging bv lEvans et al.l (i2003i) . the 10 size of the nu- 
cleus is 0'.'3 and the 25 /xm size is less than 0'.'6, which is 
the diffraction limit of the 10 m telescope. (5) The nu- 
cleus is also compact (~0'.'l-0'.'4) and has a high bright- 
ness temperature (> 10^ K) at ^2 GHz (iCondon et al.l 
[19901: IKewlev et alfeoOOt IBaan fc KlockneHl2006n . Such 
a radio nucleus is either an a ctive nucleus or a com- 
pact starburst or their hybrid ([Smith et al.l ll998'). (6) 
Chandra X-ray observations "may imply the presence 
of a Compton-thick AGN" although the identification 
is "somewhat ten tative" because o f the limited pho- 
ton statistics (M aiolino et al.l I2003D . (7) NGC 4418 
has far-IR excess compared to its radio strength, at- 
tributed either to a dust-enshrouded AGN or to a young 
compact starburst ([Kawara et al.|[l990l : lYun et ahl 120011 : 
iRoussel et aIll2003D . (8) Molecular gas concentration to 
the central < 5" has been observed with interferometers 
jlmanishi et al.l l2004t iDale et al.l l2005t iSakamoto et al.l 
[2010; C12). The mass of the molecular gas will be about 
1 X 10^ if the av erage conversion factor in our Galaxy 
([Hunter et al.lll997| ) a pplies to the observed CO(l-O) flux 
of 10q±50 Jv km s~^ feawara et al.lll99tt ISanders et al 



T99l |lYoung et al.lll995HDaie' et a .1 120051 : lAlbrechT eta 



20071 ). NGC 4418 has a high infrared luminosity-to- 
molecular gas mass ratio of LTRZ-Afmni ~ 100 Lp)IMp, 
for th e above Mmoi ([Kawara et al.l [l990l : ISanders et al.l 
119911) . High-excitation molecular lines with excita- 
tion temperatures ab ove 100 K have been detected 
towa r d the nucleus ([Aalto et al.l [2 007; Lahu is et alJ 
^007; Costagliola & Aaho 2010; Sa kamoto et al.l 120101 : 
Gonzalez- Alfonso ct al. 2012a; C12). The last four pub- 
lications report vibrationally excited molecules and the 
last three report lines involving energy levels as high 
as ^1000 K. (9) Redshifted abs orption lines have been 
found toward the nucleus by iGonzalez- Alfonso et al.l 
(l2012al) in the mid-infr ared O I and OH lines and by 
[Costagliola et al] ([20121 ) in the H 1 21 cm line. These are 
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attributed to a gas inflow. 

2. SUBMILLIMETER ARRAY OBSERVATIONS 

2.1. 350 GHz (860 fim) Observations 

We observed NGC 4418 in 2009, 2010, and 2012 
at around 350 GHz (860 /im) using the Submillime- 
ter Array (SMAj^- We used the subcompact (SC), ex- 
tended (EXT), and very extended (VEX) array configu- 
rations. Observational parameters are in Table [21 The 
data from the SC conf iguration were already reported in 
ISakamoto et all (|2010l hereafter SIO) and are combined 
here with the new data. The overall coverage of baseline 
lengths is 6-509 m or 7-590 kA at 350 GHz for the galaxy. 
We simuhaneously observed CO(3-2) and HCO+(4-3) 
lines in 2009 and CO(3-2) and HCN(4-3) in 2010 us- 
ing 2 GHz of bandwidth. We used in the EXT and SC 
configurations 4 GHz bandwidth to simultaneously ob- 
serve CO (3-2), HCO+(4-3), and HCN(4-3). The SMA 
primary beam covers about 34" (5.6 kpc) in full width 
at half maximum (FWHM) at our observing frequencies. 
The SMA spectrocorrelator was configured to the reso- 
lution of 3.25 MHz (about 3 km s~^) and recorded data 
from both upper and lower sidebands (USB and LSB). 
We used 3C273 for our gain, passband, and pointing cal- 
ibrations. We also observed the quasar J1222-I-042 and 
Titan in our 2010 VEX observations for further calibra- 
tion and its verification. Flux calibration was made us- 
ing Titan and Ganymede in 2009, and Titan in 2010 
and 2012. The angular distances from our gain cali- 
brator 3C273 to NGC 4418, J1222-h042, and Titan in 
our 2010 observations, were only 3.0, 2.7, and 4.0 de- 
grees, respectively. This mitigated calibration errors due 
to direction-dependent phase de-correlation, atmospheric 
transmission, and pointing offsets. 

We reduced our data using MIR for calibration 
(IScoyille et a ll [1991 . MIRIAD for i maging (jSault et al.l 
I1995D . and AIPS teri dle &: GreisenI [l99l and stand- 
alone programs for further data analyses. The data were 
first imaged at 30 km s~^ resolution for the entire band- 
widths to identify lines. Then the 'line-free' channels 
were identified and integrated to create continuum data. 
That continuum is then subtracted from the original to 
make continuum-subtracted visibilities. The spectra and 
the identified lines are presented in §3.11 The line data 
were further analyzed after binning to 10-30 km s~^ res- 
olutions. We did self-calibration of phase using the con- 
tinuum data when making images from the VEX data, 
although we did not do that when measuring the source 
position or fitting the visibilities to determine the size 
and fiux density of the continuum source. Table [3| sum- 
marizes the properties of the images and spectral data 
cubes used in this paper. 

2.2. 660 GHz (450 ^m) Observations 

We also observed the galaxy at 660 GHz (450 /im) dur- 
ing our VEX 350 GHz observations in 2009 and 2010 by 
allocating 2 GHz of correlator bandwidth to our high 
frequency receivers. Only the 2010 run resulted in use- 
able data from five antennas. Atmospheric absorption 

® The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics, and is funded by the Smith- 
sonian Institution and the Academia Sinica. 



was significant at our observing frequency even though 
the weather condition was excellent for Mauna Kea. The 
mean zenith opacity was 0.038 at 225 GHz and 0.89 at 
350 /im during our observations according to the mon- 
itoring at the adjacent Caltech Submillimeter Observa- 
tory. The opacity should be about the same at 450 and 
350 fj,m. The median double sideband (DSB) system 
temperature was about 1200 K toward NGC 4418. We 
used 3C273 for both passband and gain calibration. Both 
350 GHz and 660 GHz data were used for the gain cal- 
ibration. We first applied the 350 GHz gain solution 
with its phase multiplied by the ratio of the observing 
frequencies, 1.9, and then derived another gain solution 
using the 660 GHz data and applied it to improve the cal- 
ibration. No strong line is in our passband and the data 
are treated as continuum. Since Titan is too heavily re- 
solved, flux calibration at 660 GHz is based on the 350 
GHz flux density of 3C273 scale d to 660 GHz using the 
spectral index of a = —1.0 (ISlOf ). We estimate our flux 
measurements at 660 GHz to be accurate only to 36% 
(for a{a) = 0.3 and 30% error for other calibrations). 

3. OBSERVATIONAL RESULTS AND ANALYSIS 
3.1. Spectra and Line- Continuum Separation 

Figures [2l and [3l show spectrgQ made from the VEX 
data alone and from the EXT and SC data, respectively. 
The SC-I-EXT data are two times more sensitive than 
the SC data alone reported in lSlOl 

Two lines among about ten in the spectra are 
new detections in this galaxy. Both are from vibra- 
tionally excited HC3N. Namely, HC3N(z;7=li-'',J=39- 
38) is at /rost=356.072 GHz and has its upper energy 
level at = 663 K. HC3N(^;6 = ll^J=39-38) is at 
/,^,t=355.278 GHz and has = 1059 K. The dou- 
blet counterpart of the former line (117=1 ^'^,J=3 9-38) 
was already identified at /rest=355.566 GHz in ISIOI . 
Lo wer J transitions o f the l atter line were detected 
in iCostagliola fc Aahol (|2010f ) and iCia Its doublet 
counterpart (?;6=l^^,J=39-38) at /rest=355.557 GHz 
is blended with the HC3N(i;7=li^J=39-38) mentioned 
above. Although unidentified, there may be more lines, 
e.g., at /rost ~ 343.4 and 357.5 GHz; species having a 
line at the former frequency include H2CS and HOCO+. 

The 350 GHz continuum in each sideband was 
made by averaging the following channels. For the 
VEX data we used channels outside of 1900-2300 
km s"i for CO(3-2), HCN(4-3), H"CN(4-3), HC15N(4- 
3), HCN(w2=li-^,J=4-3), HCO+(4-3), CS(7-6), and 
HC3N(39-38) and (38-37). For the EXT and SC data we 
excluded channels having velocities of 1850-2350 km 
for the same linefi plus HC3N(w7=li-'^,J=39-38) and 
(w7=l^'^,J=39-38), channels between CO(3-2) and its 
nearest band edge, and the same for HCN(4-3). The 
660 GHz continuum was averaged over the full band- 
width and both sidebands. 

^ Throughout this paper, velocities are with respect to the local 
standard of rest (LSR) and are defined with the radio convention, 
i.e., Dradio = c(l — (/oba//rost )) where / denotes frequency and cthe 
speed of light. The conversion to velocity in the optical definition 
is fopt = cz = riradio + 15 km s~^ for NGC 4418. 

* The velocity range used here is the same as the one used 
to analyze the SC data in lSlOl although the note in their Table 1 
described it otherwise. 
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3.2. Concentration of Gas and Dust in the Nucleus 

Submillimeter emission is detected only from the cen- 
tral kpc (~6") with different spatial extents in various 
emission components. We show this below using maps, 
aperture photometry, and visibility fitting. 

3.2.1. Maps 

Figure |4] shows CO (3-2) channel maps to display faint 
extended emission at 0'.'6 resolution. The data cube 
was made by combining all our data. We detected here 
797±80 Jy km s^^ or 84±18 % of the single-dish C0(3- 
2) flux observed in the central 14" of the galaxy by 
lYao et all (|2003[ ). Despite the high recovery rate, the 
CO (3-2) emission was detected only in the central 6" (1 
kpc) of the galaxy, showing that the extent of the C0(3- 

2) emission is intrinsically smaller than our 34" (5.6 kpc) 
field of view. The CO emission is not only compact but 
also strongly peaked toward the nucleus within the cen- 
tral kpc. In this data set more than half (58 %) of the 
CO (3-2) flux in the central kpc is in the central 165 pc 
(1", FWHM). 

Figures E] and E] respectively are CO(3-2) and HCN(4- 

3) channel maps at ~0"3 (50 pc) resolutions showing 
that the line emission is strongly peaked at the nucleus 
at subarcsecond (<100 pc) scales. These data cubes were 
made from the VEX data alone. This CO (3-2) data cube 
contains about 50% of the CO flux detected in our 0"6 
resolution data. 

Our continuum and line moment maps are shown in 
Figs. [7] and [5] The spatial resolutions are about 0'.'6 and 
0"3 for 350 GHz data and 0"2 for 660 GHz continuum. 
In particular. Fig. |8] (k) and (1) show our detection of a 
compact continuum source at the galactic center. The 
continuum position is listed in Table [T] and agrees to 
~0'.'3 with the optical and radio positions of the nucleus; 
see the caption of T able [T] for the optical position and 
iCondon et al.l (|1990[ ) for the radio position. 

3.2.2. Aperture Photometry 

The degree of central concentration, or compactness, 
is compared among various lines and continuum at 350 
GHz in Table|4]using the ratio of flux (or flux density) be- 
tween the central 0'.'5 (80 pc) and 6" (1 kpc). The former 
photometry is made using our data fr om t he VEX con- 
figuration and the latter data are from lSlOl that used the 
SC configuration. The ratios for various emission compo- 
nents indicate their relative compactness. The relatively 
extended nature of CO(3-2) emission is evident from its 
low ratio. Only one third of the CO(3-2) emission is 
from the central 0"5. Lines other than CO have higher 
flux ratios, and hence are more concentrated toward the 
nucleus. The ratio for 860 ^m continuum, 0.95 ± 0.07, 
is the highest among the ratios with S/N > 3 in Table 
131 Thus 860 /im continuum emission is more compact 
than both the CO and non-CO high excitation lines. At 
the same time, the continuum emission in the VEX data 
is slightly more extended than the 0"3 synthesized beam 
as seen in the last two panels of Fig. [5l The detected 
450 yum continuum also looks compact but it has only a 
minor fraction of single-dish measurement. 



Figure ini (a)-(c) show fits to the continuum visibilities 
of NGC 4418, a test point source J1222-I-042, and Titan 
whose subarcsecond size is known. The fit results are 
hsted in Table El For NGC 4418, data from the three 
array configurations are combined. This cannot be done 
for the calibrators because their fiux densities and the 
apparent diameter of Titan change with time. Since the 
NGC 4418 data were taken in three different nights the 
three data sets should have independent flux calibration 
errors, which are typically about 10-20% for the SMA 
350 GHz observations. However, we did not add the flux 
scaling error in the plot and in our fitting because do- 
ing so would make the reduced well below unity. It 
is most likely that we achieved better flux calibration 
than usual, presumably because we had accurate point- 
ing from repeated pointing on the bright quasar 3C273 
next to NGC 4418, our primary flux calibrator Titan was 
also close to 3C273 and NGC 4418, and atmospheric ab- 
sorption was low (225 GHz zenith opacity < 0.05 for all 
the three nights). The size of 850 /im continuum emission 
of NGC 4418 obtained from our data is 0"10 (FWHM). 
We will use this for the source modeling in §3.31 

We verified our visibility fitting using the quasar 
J1222-f042 and Titan. They were gain calibrated us- 
ing 3C273 in the same way as NGC 4418 was calibrated. 
Our fit of J1222+042 using a Gaussian shows that blur- 
ring of the point source was less than 50 milli arcsec 
(mas) in our observations. The diameter of Titan at 860 
is known to be 5230 km or 0'.'8425 for the distance of 
8.5588 AU. This size includes the 40 km altitude of the 
tropopause, around which is the sub-mm photosphere of 
Titan. The distance at the time of our 2010 VEX obser- 
vations is from JPL-Horizons. Fig. |9] (c) shows that our 
data match the expected visibility curve very well with 
only flux scaling. When the data were fit for both the 
source size and amplitude scaling, we obtained the size 
of the sub-arcsecond source with an accuracy better than 
3% or to 20 mas. 

The flux density of Titan sets the absolute amplitude 
scaling of all the visibility fltting. We adopt the con- 
tinuum brightness temperature of 77 K as the Planck 
temperature that gives the total flux density of Titan for 
the solid-surface diameter (5150 km) of the satellite (M. 
Gurwell, private communication) . Our Titan data do not 
contain CO(3-2) and HCN(4-3) lines, both of which are 
known to be bright. Still, our imperfect knowledge of the 
atmosphere of Titan may introduce a common amplitude 
error to our data even though the relative calibration us- 
ing Titan is precise and consistent among the three data 
sets. We therefore assign 5% to the overall flux scal- 
ing uncertainty of our NGC 4418 flt, even though the 
formal error in the fitting was about 1%. In addition 
to this, there is likely contamination of line emission to 
the 'line-free' channels used for the USB continuum of 
the galaxy. We observed a sign of that in our SC con- 
figuration data in the form of an uncomfortably large 
spectral index a = 6.25 ± 0.84 {Si, oc i/") between the 
two sidebands (ISlOl ). A similar sign was seen in our 2010 
VEX data too. Since we used the USB continuum in the 
multi-configuration fitting as LSB data were noisier, the 
line contamination is likely in the fitted amplitude. Its 
effect on the luminosity estimate is addressed in >j3.3.3l 



3.2.3. Visibility Fitting : 860 yim Continuum 



3.2.4. Visibility Fitting : Lines 
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Fig. inialso shows visibility amplitudes of line emission 
against baseline lengths. Line visibilities were averaged 
over 1900-2300 km s~^ after continuum subtraction, and 
then binned according to the baseline length. It is evi- 
dent in the figure that the line visibility amplitudes de- 
cline faster than the continuum amplitudes. In other 
words, line emission is more extended than continuum. 
Comparison of CO, HCN, and HCO+ data tells us that 
the CO emission is more extended than the latter two, 
confirming our observation using the aperture photome- 
try. Moreover, the CO visibility plot shows that the data 
cannot be fit with a single Gaussian in the u-v domain. 
Hence the spatial distribution of the CO(3-2) emission 
is not a Gaussian. The line emission distribution has at 
least a compact subarcsecond peak toward the nucleus 
(Fig. IS]) and an extended component with the total ex- 
tent of ^5" (Fig. [7]). The extended component is non- 
axisymmetric around the nucleus (Fig. [7])- 

3.3. Parameters of the Continuum Core 
3.3.1. Spectral Energy Distribution, Continuum Opacity 

The continuum SED has a power-law slope of a = 
2.55±0.18 {S^ oc z/") between 1.3 mm and 0.85 mm. This 
is measured fro m the SMA data in Figure [TUl namely 
from this work, ISIOI and the 1.3 an d 1.1 mm obser- 
vations used in the companion paper (|C12D . Since the 
SMA continuum data were taken from line- free spectrom- 
eter channels in each dataset, none of the SMA data 
is contaminated by strong emission lines such as CO 
unlike wide-band bolometer data. The spectral slope 
that we obtained is shallower (i.e., a is smaller) than 
those of other local IR luminous galaxies but is con- 
sistent wi th the r esults of the FIR-to-mm SE D surv eys 
of LiscnfcldiFil] pOOO) and Yang & PhillipsI (pOO?! ), in 
each of which NGC 4418 has the shallowest spectral slope 
in about 15 galaxies. 

Two simple models for the spectral index suggest that 
the a measured at the nucleus is due to the dust opac- 
ity T ^ 1 at A ~ 1 mm. In the first model of a sim- 
ple one-zone slab, the spectral slope should be a = 
2 -I- /3T(e'^ — 1)~^ where (3 is the index of power-law fre- 
quency dependence of the dust opacity and r is the slab 
opacity at the wavelength of the a measurement. Here 
we used the Rayleigh- Jeans approximation because the 
brightness temperature of the O'.'l nucleus exceeds 100 K 
as we will see next (i j3.3.3p . With optically thin emis- 
sion, this model gives the familiar form of a = 2 -I- /3. 
The low a observed in NGC 4418 means low /3 or a 
moderate opacity of dust emission, i.e., r ^ 1 at A ^ 1 
mm. The dust op acity a bout 1 is favored, as was con- 
cluded by Liscnfcl d et al.l ([2000), because most infrared 
luminous galaxies hav e /3 ~ 1.5-2 (jLisenfeld et al.l 120001 : 
IDunne fc EalesI [2001D . The observed a corresponds to 
T = 2.2 ± 0.5 and 1.8 ± 0.5 for /3 = 2 and 1.5, respec- 
tively. For the second model, we note that the SED slope 
does not become 2 at shorter wavelengths as expected 
for T 3> 1. The shallower slope of the Planck function 
compared to the Rayleigh- Jeans formula must partly be 
the reason. Another likely reason is that the dust dis- 
tribution is not just a O'.'l slab of uniform opacity but 
presumably has a low opacity halo as suggested by the 
larger extent of line emission than the 860 fim contin- 
uum. The spectral slope in the general multi-component 



case without mutual shielding is the flux-weighted mean 
of the spectral indexes, a — J2 Siai/Y^ Si, where the ith. 
component has the fiux density Si and spectral index a^. 
For ai — 2 and a2 = 4, the observed a corresponds to a 
fraction of the optically thick core of 73 ± 9 %. Thus as 
long as we assume /3 w 2, we need an optically thick core 
(r > 1) that dominates the spectral energy distribution 
at A ~ 1 mm. In addition to the two models, the high 
opacity is also supported by the high brightness tem- 
perature that sets a lower limit of about 0.1 to the 860 
fim opacity from the constraint that opacity-corrected 
brightness temperature cannot exceed dust sublimation 
temperature. 

There are other factors that depend on frequency and 
potentially affect the spectral index. They include the 
core size defined with the t — 1 surface, the dust tem- 
perature at the photosphere, and absorption of the core 
emission by dust around it. The last one is addressed in 

3.3.2. Gas Column Density, Mean Density, and Mass 

The dust continuum opacity Tseo /xm ~ 1 of the core 
translates to a hydrogen column density Nn ^ 10^^'^ 
cm~^ for the dust opacity to column density relation 
N(R + }i2)/Tx = 12 X 10^^ (A/400 ^m)^ H-atom cm'^ 
(jKeene et al.lll983lHildebran d 1983). Since the core has 
a size of ~20 pc (O'.'l in FWHM from > j3.2.4l see the next 
subsection for two models), we obtain riHs ~ 5 x 10^ 
cm~^ for the mean gas density (i.e., number density of 
hydrogen molecules) and M„^oi{r < 10 pc) ~ 1 x 10* Mq 
for the total gas mass in the core assuming that the gas is 
mostly molecular there. The total mass depends on the 
density distribution. The estimate above is for a uniform 
density and is an upper limit if the gas density increases 
toward the center. 

3.3.3. Brightness Temperature and Luminosity 

Table [6] lists the deconvolved (peak) brightness tem- 
perature and the bolometric luminosity of the nucleus. 
The luminosity is derived from the 860 fim core size and 
the temperature. Here we use two different models to 
estimate the brightness temperature a nd th e bolometric 
luminosity of the nucleus as we did in IS08I for Arp 220. 
One is the circular Gaussian we used for the visibility 
fit above and the other is a circular disk with uniform 
brightness. The Gaussian model describes dust that is 
not fully opaque across the nucleus and that may also 
have a radial temperature gradient. The circular disk 
model corresponds to a sphere of uniform surface bright- 
ness projected on the sky. Both models can fit the same 
data well when the visibilities are sampled only within 
the half-maximum baseline length of the visibility func- 
tion. The diameter of the model disk is 1 .6 times the 
FWHM of the Gaussian in such a case ([SOl). For the lu- 
minosity calculation, it is assumed that the 860 /im emis- 
sion is optically thick in the uniform disk/sphere model. 
Since the deconvolved Tb was used without any opacity 
correction (i.e., not divided by 1 — e""^), the derived tem- 
perature and luminosity are lower limits. We calculated 
the luminosity of the Gaussian model by using its FWHM 
size and its peak brightness temperature. This crudely 
approximates the effect mentioned earlier that the size of 
the photosphere should increase toward higher frequen- 
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cies as dust opacity increases with frequency. 

The derived luminosity of the nucleus is about 1 x 
10^^ Lq, which is comparable to the total luminosity of 
the galaxy (10^^ "'^ L©, Table [T]). We obtained nuclear 
luminosities log(iboi/-^o) — 11.3 ±0.3 and 10.8 ±0.3 for 
the two models. Luminosity surface density and volume 
density are also calculated for each model and listed in 
Table [6] for the central ~20 pc. The errors here and in 
Table [6] are random errors due to the 5% flux error and 
the 11% size uncertainty and do not include systematic 
error due to the models. 

There are two notable sources of systematic errors in 
the derived parameters in addition to the assumptions on 
the source shape and on the 860 /im opacity. (The under- 
estimating effect of the latter is minor as noted in SOS.) 
One is the absorption of the core emission by the sur- 
rounding dust and gas. It arises because the ^O'.'l ('--^20 
pc) continuum core is embedded in an order of magnitude 
larger envelope as we saw in tj3.2l The continuum emis- 
sion from the core must suffer from extinction to some ex- 
tent by the envelope through continuum self-absorption 
and line absorption by various molecules. Our continuum 
brightness temperature derived without the extinction 
correction is therefore a lower limit of the true bright- 
ness temperature at the surface of the ~0'.'l core. The 
correction for it would increase the core luminosity. We 
will analyze this submillimeter extinction a little more 
in §3.3.41 and find the extinction to be about 5% at 860 
/xm. The other source of systematic error is the possible 
contamination of the 860 /im continuum by molecular 
lines. This effect can be estimated to be about 10% as- 
suming that the spectral index between lower and upper 
sidebands obtained in the SC configuration, 6.25 ± 0.84 
( §3.2.3|) . is due to line contamination to the USB data 
and that the true continuum spectral index is 2.55 ±0.18 
( §3.3.1|) . The continuum brightness temperature of the 
core would be lower without the contamination. How- 
ever, this does not necessarily mean that the bolomet- 
ric luminosity of the core is overestimated by this effect. 
Thermal lines can significantly add flux density to the 
860 /im continuum only when the line and continuum 
opacities satisfy runc > Tcont and Tcont ^ 1, because no 
thermal line will be seen from a blackbody. Since we use 
the brightness temperature of the continuum core for our 
luminosity estimate assuming that the continuum is op- 
tically thick, the line emission, if it is from the O'.'l core 
itself, would (partly) compensate for the decrement of 
the observed continuum brightness temperature due to 
Tcont that is not ^ 1. 

In addition to these two error sources, our visibility 
fits are insensitive to a point-like source with a small 
flux contribution. This is due to our u-v coverage limit. 
For example, a 20 mJy source with a size of 20 mas (3 
pc) can be undetected in our visibility flt (Fig. [9] a) and 
would have a temperature of 700 K and luminosity of 
10"'^^ Lq. The total luminosity of the galaxy limits the 
presence of such compact sources. 

3.3.4. ^50 nm Continuum 

We detected 450 /im continuum at 5.8a with the res- 
olution of 0'.'23 X 0'.'15 at the position of the 860 /Ltm 
continuum peak (Fig. [5] (1)). The 450 fj,m. nucleus in 
the map appears unresolved but this is largely due to 
the low signal-to- noise (S/N) ratio. (Image-domain fit- 



ting gives a nominal deconvolved size of 0'.'10±0'.'04 in 
the major axis while the source is unresolved in the mi- 
nor axis.) The 450 /im detection of the subarcsec-scale 
nucleus is at lOu and hence even more secure when the 
data are convolved to 0'.'5 resolution (Table SJ. (The 
S/N increases because the convolution lowers the weight 
of noisy long baselines.) The peak and total flux den- 
sities in both resolutions are about 370 mJy and mJy 
beam~^, respectively. 

The 450 /im continuum flux density that we detected 
is signiflcantly lower than that of the central 3 kpc of 
the galaxy (see Fig. [IHl). The ratio is 0.31 ± 0.15 
between the SMA and single-dish mea s ureme nts. The 
latter was made by iRoche fc Chandlell ()1993l ) with the 
UKT14 bolometer, the filter of which had the effective 
frequency of 682 GHz (440 um) and the bandwidth of 
84 GHz (jPuncan et al.l il99oV We subtracted from the 
bolometer measurement the contribution of CO(6-5) line 
(~0.1 Jy, c alcula ted from the CO (3-2) single-dish flux of 
lYao et al] ()2003[ ) scaled by a factor of 4) and rescaled the 
continuum to our observing frequency using the spectral 
index a = 2.55. The single-dish flux density of the 450 
fim continuum estimated this way is 1.18 ± 0.33 Jy. The 
low flux recovery ratio at 450 /im suggests that there is 
a significant 450 /xm emission that is too extended for 
detection at our sensitivity and with the shortest base- 
line of 240 kA in the VEX configuration. We simulated 
the degree of flux recovery with our actual u-v coverage 
for Gaussian sources of various sizes. (We did not fit 
the visibilities unlike our 860 /im data analysis because 
the 450 fim data have lower S/N ratio.) Fig. [TT] shows 
that the size of the 450 /im continuum source would be 
0'.'24 ± 0'.'07 (FWHM) if the source were approximately a 
single Gaussian. The peak brightness temperature of the 
source would be 70 ± 20 K and the bolometric luminosity 
of this source would be on the order of only 10^°'^ Lq, 
which disagrees with our estimate at 860 /im. 

A two-component model with a core and an envelope 
agrees with the observations within the errors and is 
more reasonable in our view. The two components are 
introduced in view of the compact continuum core and 
the more extended CO emission both observed at 860 
fim. They are treated separately to simplify our model- 
ing though they are probably not two discrete entities. 
( More detailed multi-comp o nent m odel s have been made 
bv lGonzalez- Alfonso etldl (|2012af ) and lC12l Our inten- 
tion here is not to update their models but to apply a 
simple model to understand our sub-mm continuum ob- 
servations.) The core and envelope are assumed to have 
sizes of O'.'l and 0'.'5 (both in FWHM), respectively. Also 
the core is assumed to dominate the 860 /im emission 
with 85% of the total fiux density. The extended enve- 
lope with the small remaining flux density has too low a 
surface brightness to be detected in our high-resolution 
map. However, the amount of extended 860 /im con- 
tinuum in this model is implied by our CO (3-2) data. 
Namely, if one third of the CO (3-2) emission is from the 
envelope region, then its fiux of about 300 Jy km s~^ and 
the typical CO(3-2) equivalent width of 1 x 10"^ km s~^ 
among infrared bright galaxies (jSeaauist et al.|[20M) sug- 
gest a 860 /im extended continuum of ~30 mJy, which 
is about 15% of the total co ntinuu m flux density. The 
analysis of spectral index in i j3.3.1l also implied a small 
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contribution to 860 /zm emission from an optically-thin 
and extended component. We assume the spectral index 
between 860 fj,m and 450 /im to be 2 for the optically 
thick core and 4 for the optically thin envelope (i.e., the 
model assumes /3 — 2). In addition, since the core is 
embedded in the envelope, the former should be extin- 
guished by the latter. The opacity of the envelope to 
cover the core is assume d to be ^1 at 200 iim, which we 
took from the model of iGonzalez- Alfonso et all ()2012af) 
noting that it makes the envelope optically thin at 860 
/im (r — 0.05) as we assumed above. At 450 fim the core 
is extinguished by the envelope with r = 0.2 and the en- 
velope itself is resolved out in our VEX data. Overall, 
this model predicts our flux recovery rate to be about 
45%, which agrees with our observed value of 31 ± 15 %. 
The current data quality does not justify further param- 
eter tuning and model refinement. 

This model illustrates that two effects likely contribute 
to the lower flux recovery rate at 450 /im than 860 /im. 
One is that the core emission is saturated while emission 
from the extended outer region is not. The other is that 
the emission from the hot nucleus is extinguished more 
at shorter wavelengths by the colder circumnuclear dust. 
These effects can also explain the 60 /im peak of the 
SED of NGC 4418 despite the energetically dominant 
core emitting at > 100 K at 860 /im. The same effects 
probably explain similar observations in Arp 220 where 
the two nuclei contr ibute much less to the total at 435 
fim than at 860 fj,m ([Matsushita et aLll2009f ). 

3.4. Kinematics of Molecular Gas 

The HCN velocity map Fig. ^i) shows in the central 
0'.'5 (80 pc) a systematic gradient of mean velocity in 
the northeast-southwest direction, which is close to the 
galaxy major axis (p. a. =60°, Table[T]). The HCN channel 
maps in Fig. |6] show this velocity structure as the shift 
of the emission peak across the continuum position from 
southwest (lower-right) to northeast (upper-left) in chan- 
nels between 1990 and 2200 km s'^ In Fi gure 1121 mo- 
ment maps of HCN, HCO''", and CS confirm the presence 
of the gradient across the nucleus although the direction 
of the velocity gradient is closer to the north-south di- 
rection in these lower resolution data. This is also the 
case in the CO high-resolution data (Fig. [8] (b)) in the 
central 0'.'5. 

We attribute the central velocity gradient to rotation 
around the luminous core and estimate its dynamical 
mass to be M^y^^r < 15 pc) ~ 2 x IO^Mq. We estimated 
for this a velocity shift of about 400 km s~^ in 0'.'2 across 
the nucleus from the CO and HCN position- velocity dia- 
grams in Fig. 1131 where almost the full velocity range of 
each line is observed at the central ~0'.'3. The Keplerian 
dynamical mass above is for the 62° inclination of the 
galaxy. Because of the insufficient spatial resolution to 
obtain an accurate rotation curve, our Mdyn estimate is 
admittedly crude and its error may be as large as a factor 
of 2 even if the adopted inclination is correct. If the rota- 
tion velocity increases toward the center at these scales, 
as we will see possible below, the dynamical mass is more 
likely to be overestimated. For comparison, we estimated 
in our companion paper the mass of the central black 
hole in this galaxy to be 9 x 10^ Mq (jCostagliola et all 
120121 for the distance of 34 Mpc). This value follows 



from the K band bulge luminosity 10^'^ Lq^k, which is 
calculated from the total galaxy luminosity in Table [T] 
and the me an bulge-to-total luminos ity ratio of 0.31 for 
Sa galaxies (iGraham fc Worlevl[200l . The mass should 
have a factor of 2 uncert ainty due to the scatter in the 
Lk, huigc-Mi,i, correlation ([Marconi fc Huntll2003D . 

We also infer from information of multiple lines that 
gas motion is faster toward the center. Fig. [T4]and Table 
[7]show the centroid velocities and widths of the molec ular 
lines measured in the SC configuration data of 'SlO. An 
interesting observation here is that there are two groups 
in the centroid velocity - line FWHM plane. The first 
group consists of CO(3-2), C0(2-l), and N2H+(3-2) and 
have centroid velocities Vc ~ 2100 km s~^ and FWHM w 
142 km s~^. The rest of the lines form the second group 
at Vc « 2088 km a'^ and FWHM w 234 km s'^ Also, 
there is a trend that lines with higher critical densities 
are wider. This variation of line parameters again sug- 
gests that the molecular gas in the center of NGC 4418 
has non-uniform properties or multiple components. In- 
deed, the latter trend can be easily explained if the gas 
motion around the nucleus is faster closer to the center 
and if denser gas is more localized around the nucleus. 
This model agrees with our previous observation in Ta- 
ble H] that the lines with higher critical densities are more 
concentrated in the nucleus than CO. If lines with higher 
critical densities indeed trace gas closer to the nucleus, 
and if we also assume that the gas distribution is axisym- 
metric around the nucleus, the true systemic velocity of 
the nucleus would be about 2088 km s~^. 

We can also evaluate radial motion of molecular gas 
from our data. HCN spectra near the continuum peak in 
Fig. [15] and the HCN PV diagram (Fig.[T3|) show a minor 
depression {2a) at 2080 km s~^. This may be due to gas 
at about Vgys and through line self-absorption or line ab- 
sorption of the continuum from the core. The velocity of 
this possible absorption and the absence of clear absorp- 
tion off the systemic velocity indicate little radial motion 
of the HCN gas. Meanwhile, the CO emission is signifi- 
cantly brighter below about 2100 km s~^ than at higher 
velocities in the PV diagram. See also the CO line profile 
at the continuum core position in Fig. [15] where the line 
profile is asymmetric. If this is due to self- absorption 
of the CO emission, the foreground gas causing the ab- 
sorption (and weaker CO emission) is redshifted with 
respect to the systemic velocity and hence is flowing in- 
ward. Although this asymmetry can be due to a chance 
imbalance of gas distribution between approaching and 
receding halves around the galactic center, this is quali- 
tatively consistent with the redshifted O I, OH, and H I 
absorption obs erved bv IGonzalez- Alfonso et all ([2012af) 
and Costagli ola et al.l ([2012t ). The plausible gas motion 
mentioned here is opposite to that in Arp 220 where P- 
Cygni pro files of molecular li nes suggested molecular out- 
flows (Sak amoto et al.ll2009D . 

Gas motion outside the central 0'.'5 does not appear to 
be ordinary rotation in our data. The CO mean velocity 
maps Figs.[7|(b) and[8](b) do not show an ordinary spider 
pattern of a rotating gas disk at r > 0'.'25. The velocity 
maps probably contain a mixture of rotational and radial 
motions and appear complex because of non-uniform gas 
distribution and optical thickness of the CO line. 
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3.5. Gas Properties 

Our high resohition observations also provide the fol- 
lowing new insights into gas properties in the nucleus. 

First, the CO(3-2) line has a peak brightness tempera- 
ture of 90 K at (X'S resolution (Tabled]), indicating warm 
molecular gas around the nucleus. The brightness tem- 
perature of thermalized CO is the gas kinetic temper- 
ature diluted by the CO opacity and by the observing 
beam if it is larger than the source. Both dilution effects 
are probably small for the extended ^^CO in this case. 
To put the observed Tb into context, the peak value of 90 
K at 50 pc resolution surpasses the peak CO brightness 
temperatures of about 50 K observed in Arp 220 and 
NGC 253 at 100 pc (O'.'S) and 20 pc (1 ") resolutions, re- 
spectively (Safemoto_eiaiJEQOllMll!)- Meanwhile both 
HCN(4-3) and HCO+(4-3) have peak brightness temper- 
atures of about 30 K at the same 0'.'3 resolution (Table 
|4]). Their ~3 times lower Tb than that of CO must be 
largely due to beam dilution effects for these lines be- 
cause they were found more compact than CO (3-2) in 
i j3.2l Subthermal excitation and lower opacities of HCN 
and HCO^' may be additional reasons for the lower Tb 
of these lines. 

Second, the detection of lines from vibrationally ex- 
cited molecules is another indication of high tempera- 
tures in the nucleus. The HCN and HC3N lines that 
we detected from their vibrationally excited states have 
upper energy levels 510-1070 K above the ground level. 
Vibrational temperatures of these molecules have been 
estimated to be around 300 K ([Costagliola fc Aaltoll2010l : 
iSlOtlCl^ . Infrared radiation has been suggested to play 
a significant role for the vibrational excitation in the 
~100-200 K gas and dust that we inferre d fro m bright- 
ness temperatures. Our companion paper ()C12f ) presents 
a more detailed analysis and discussion on the excitation. 

Third, the mean gas density tt-Hs ~ 5 x 10^ cm"'^ of 
the ~20 pc core (i j3.3.2l) is high enough for CO excita- 
tion to J=3 but not for HCN and HCO+ to J=4. Their 
critical densities for coUisional excitation at 100 K are 
~10'*'^, 10^, and 10^ cm~^, respectively. The mean den- 
sity is also far short of the '^lO^^ cm~^ critical density 
to vibrationally excite HCN with H2 collisions. There- 
fore the excitation of the lines with high critical den- 
sities are probably achieved through photon trapping, 
non-uniform density distribution to allow higher local 
densities than the meari, mid -infrared radiation, and/or 
electron colhsions (|S10l : lC12f ). 

Fourth, the ratios between lines and between line and 
continuum vary with the area size of the measurements in 
the way expected when the interstellar medium is warmer 
and denser at smaller radii. This is a consequence of 
the size variation among various emission lines fi i3.2p . 
For example the HCN(4-3)/CO(3-2) ratio of integrated 
brightness temperatures is 0.32±0.03 and 0.170±0.004 at 
0'.'5 and 6" apertures, respectively. The ratio is expected 
to be higher when molecular gas is denser, warmer, or fa- 
vorable for HCN line emission in other ways. Similarly, 
the equivalent widths of lines with respect to the contin- 
uum depend on the size scale. The CO(3-2) equivalent 
width changes by a factor of 3, between 1.4 x 10'^ and 
4.2 X 10^ km s-\ from 0'.'5 to 6" scales. The smaller 
equivalent width at the smaller scale can be explained 
by the higher opacity of the dust continuum at the con- 



tinuum core ()S08[ ). Finally, the ratio between HCN(4-3) 
to HCO+(4-3) is 2.08 ± 0.27 in the central 0'.'5 while it 
is 1.65 ± 0.07 in the central 6" (fsToh . Among other dif- 
ferences, HCN has an order of magnitude higher critical 
density for collisional excitati on a nd is also easier to ra- 
diatively excite than HCO+ (|sTof ). See gXUfor more 
discussion on this line ratio. 

3.6. Optical Color Index 

We found a U-shaped red feature in the multi-band 
optical images of the galaxy from the Sloan Digital Sky 
Survey (SDSS). Fig. [TOl shows the distribution of a flux 
ratio involving four SDSS bands, {i' + z') / {g' -\-r'). Areas 
brighter at longer wavelengths (i' and z') are shown in 
red while areas of lower flux ratios are in blue. The U- 
shaped red feature is along the northwestern semi-minor 
axis of the galaxy. No counterpart is visible around the 
southeastern semi-minor axis. The reddest point is al- 
most on our 860 /Lim continuum peak (i.e., the nucleus) 
and is about 0.5 mag (i.e., a factor ~1.6) redder than 
the off-center region outside the U-shape. This feature 
extends at least up to 10" (1.7 kpc on the sky) from the 
nucleus. The redder part of it near the nucleus has an 
elongation along the position angle of about 15°, similar 
to the faint CO (3-2) emission in Fig. [71 

We interpret the red feature to be an outflow cone em- 
anating from the galactic center along the rotation axis 
of the galaxy. The feature is very unlikely to be a struc- 
ture on the galactic plane because such a feature would 
be torn by differential rotation and also because the fea- 
ture is on the minor axis of the galaxy. The minor axis 
is determined by our viewpoint while a radial feature in 
the galactic plane should have no preferred position an- 
gle with respect to us. The outflow cone model does not 
have this problem because any feature along the rota- 
tion axis of the galaxy is always projected onto the mi- 
nor axis. Also, the outflow can be bipolar in this model 
even though we only see half of it in the color index 
map. The fact that we see the red feature only on the 
northwestern side suggests the side to be the near side of 
the outflow and the far side of the galaxy disk, against 
which the outflow cone is visible as a red silhouette due 
to color-dependent extinction. Dust in the cone probably 
causes the extinction and the reddening. Our CO veloc- 
ity data in Fig. [7] cover only a small fraction of the red 
optical feature and do not reveal the gas motion in the 
red cone. The data only tell that the velocity held within 
a radius of ^^500 pc is more complicated than expected 
from purely circular rotat i on ar ound the center. How- 
ever, |LehnerF&3ickma^ ()1995D reported a 200 km s~^ 
velocity shift along the minor axis from their [N II] line 
measurements for a survey of starburst winds. The veloc- 
ity shift suggests radial motion of the ionized gas along 
the polar axis. Although we cannot yet conflrm an out- 
ward radial motion because the direction of the [N II] 
velocity shift was not given, the velocity shift along the 
red cone better fits the outflow interpretation. For this 
radial motion to be an inflow, gas and dust need to be 
falling in the U-shape cone toward the nucleus along the 
polar axis, which seems too much of a coincidence. Our 
outflow model predicts that the gas on the northwestern 
cone should be blueshifted on average with respect to the 
nucleus. 

Much closer to the nucleus, lEvans et al.l (|2003l) found 
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radial dark lanes out to about 3" from the nucleus in their 
near-infrared data. The lanes have an average position 
angle of about 15° and are on both sides of the nucleus. 
The fainter lane to the north is in the same direction 
as the reddest part of the optical U-shaped feature as 
well as the CO elongation. There might be therefore a 
change of the outflow direction, first at P. A. ^ 15° in the 
innermost region (r < 3" « 500 pc) and then along the 
minor axis of the galaxy (P.A. ^ —30°) outward. 

We estimate the gas mass in the outflow to be ^4 x 10^ 
Mq. This is from the reddening and assumes a bipo- 
lar structure of the outflow. The excess color index in 
g' — r' is 0.8 mag at the central peak, about 0.5 mag in 
the reddest 1 kpc^ near the base of the U-shape, and 1 
mag kpc^ when integrated over 2x2 kpc^ encompass- 
ing the entire U-shape. Adopting the Milky- Way extinc- 
tion law of Ag> - A^/ = OJlAv and I](H-hH2)/Ay = 
2 X 10^^ H cm~^ = 16 M© pc~^, we estimate the gas 
column density to be A^h = 2 x 10^^ cm~^ at the peak, 
the mean mass surface density in the reddest 1 kpc^ to 
be about 10 Mq pc~^, and the gas mass in the one side 
of the outflow (i.e., the U-shape) to be about 2 x 10^ 
Mq. The mass in the entire outflow is twice the last 
value on the assumption of a symmetric bipolar struc- 
ture. This mass estimate assumes foreground-screen ex- 
tinction. The 860 fim dust core mixed with stars in the 
galactic center does not contribute much to the redden- 
ing and, appropriately, to the mass estimated here. Note 
that the reddening data (and hence the outflow mass esti- 
mate from it) is more sensitive than our CO data in which 
the U-shaped cone is not detected. The cone would have 
a CO line brightness temperature on the order of only 0.1 
K when the outflow gas is 100% molecular, its CO-to-H2 
conversion factor is 2 x 10^° cm~^ (K km s"^)^^, and the 
line has a velocity width of 30 km at each position. 
This is about 0.5a in our 0'.'6 resolution CO (3-2) data 
cube. 

We further estimate the kinematical age of the outflow 
to be on the order 10 V2qq Myr and the gas outflow rate 
~4 V200 Mq yr~^. The parameter V200 is the outflow ve- 
locity normalized by 200 km and is about unity for 
the projected minor-axis velocity shear of 200 km s^^ 
measured in [N II]. The age (i.e., crossing time) of the 
outflow is the ratio between the de-projected size of the 
outflow cone, 2 kpc, and the outflow velocity. The mass 
outflow rate from the nucleus is from the gas mass above 
and the crossing time. We caution about the uncertain- 
ties in these parameters due to the unknown velocities of 
various outflow media, the poorly constrained extent of 
the flow, and the uncertain projection effect. For exam- 
ple, the ionized gas traced by [N II] and the dust causing 
the reddening may have different bulk velocities. Also 
the opening angle of the outflow and the possible bend 
of the outflow direction add uncertainty to our outflow 
velocity and to the parameters dependent on it. For ex- 
ample, if the outflow velocity is as la rge as 1000 km s~^ 
seen in some mo lecular outflows (e.g.. lChung et al.ll2011l : 
lSakamotoll2012t ). the crossing time will be reduced to 2 
Myr and the outflow rate will be ^^20 Mq yr~^. Another 
caution is that, with the velocity information currently at 
hand, we cannot tell whether the gas and dust blown out 
from the nucleus to > 1 kpc above the galactic plane will 
leave the galaxy or fall back to it. In the latter case, the 



mass loss rate of the galaxy is smaller than the outflow 
rate from the nucleus. 

4. DISCUSSION 
4.1. The Luminosity Source at the Nucleus 

We have obtained several parameters to constrain the 
nature of the luminosity source such as the iboi, M^ya, 
size, and A'h of the luminous core, and E(Lboi) 
10^-^ Lq pc"2 Lboi/Mdyn 500 Lq Mq"^ calcu- 
lated from them. While we derived these from our 860 
fim observatio ns alone they qualit atively agree with the 
observations of lRoche et al.l ()1986[ ) and many others that 
most of the luminosity of the galaxy must come from a 
compact and deeply enshrouded nucleus. The param- 
eter values are comparable to or higher than thos e in 
the ultraluminous infrared galaxy Arp 220, where IS08I 
obtained Lboi/A^dyn > 400 Lq Mq^^ and i;(iboi) > 
10^'^Lq pc^^ in the central 80 pc of the western nucleus. 

4.1.1. Constraint from L/M 

A young starburst can have a bolometric luminosity to 
mass ratio at or above 1000 Lq Mq~^ for 3-10 Myr de- 
pending on the initial mass function while a bare AGN at 
the Eddington luminosity has L/Mbh ~ 10^'^ Lq Mq~^ . 
Fig. [17] sh ows this for starbu r sts simulated using 
Starburst9 9 (iLeitherer et ap 119991 iVazquez fc Leithereri 
[2005; Lcit herer et al.l I2010D and the initial mass func- 
tion flMF) of iKroupal (|2002l ). The duration of L/M > 
1 X lO'^ Lq Mq~^ is about 5 Myr for the normal (i.e., un- 
truncated) Kroupa IMF in both instantaneous and con- 
tinuous starbursts. In a continuous starburst, the com- 
bination of L/M « 10^ Lq Mq~^ and L « 10^^ Lq is 
obtained with the star formation rates (SFRs) of 30-100 
Mq yi:~^ at the ages of 3-1 Myr. The starburst can be 
older and SFR lower if the IMF is biased toward high 
masses compared to the standard one (e.g., age^lO Myr 
and SFR- 10 Mq yr"^ for the IMF mass range of 1- 
100 Mq). The mass used for L/M in Fig. [13 is that of 
newborn stars (or the black hole mass for the Eddington 
limit). The dynamical mass in our observed iboi/Afdyn 
includes the gas covering the nucleus, the central black 
hole, and the old stellar population preexisting the star- 
burst or nuclear activity. The gas mass is not negligible 
since it was estimated to be about half of the dynami- 
cal mass fi )3.3.2p . Removing these from the denomina- 
tor of the observed L/M ~ 500 Lq Mq" will increase 
the ratio to >1000 Lq Mq~^. Fig. [ITj therefore suggests 
that if a starburst is responsible for the core luminos- 
ity of NGC 4418 its L/M is near the maximum that a 
starburst can have only when it is young, assuming the 
standard IMF. The L/M data from our submillimeter 
observations therefore provide a new way to constrain 
the age in the starburst model. 

Regarding whether a compact and massive star clus- 
ter mentioned above is able to form, it has been ar- 
gued that the maximum L/M for a cluster-forming 

gas cloud or disk is 50 0-1000 Lq Mq"^ (ScoviUc 200^; 
[Thompson et al.ll2005f ). This limit is because too much 
luminosity will blow out the star forming gas with radi- 
ation pressure on dust and halt the star formation. The 
observed L/M is c lose to this limit but does not clearly 
exceed it. iMurravl 12009,) predicted a mass-size relation 
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for massive clusters taking the radiation pressure into 
account. The diameters of 10^'^ and 10® Mq clusters, 
which can have ~10^"'^ Lq when young, are 5 and 10 
pc respectively in his relation. Thus at least these con- 
siderations do not reject the starburst model with the 
parameters we currently have. 

If the energy source of the IO^^Lq core is an ac- 
creting black hole then its Eddington ratio must be 
log(i/LEdd) = —0.5 ± 0.3 assuming the black hole mass 
of 10^ °±°-3Mq (hence LEdd = lO"-^*^-^ Lq) estimated 
from the bulge luminosity ()C12| . SjSJ]). Thus an AGN has 
no problem to be the dominant energy source as long as 
the L/Mbh is concerned. It has no problem to be smaller 
than the 20 pc core either. 

4.1.2. Constraint from E(Lboi) 

The luminosity surface density that we obtained, 
I](Lboi) = 10®-5=t°-5 Lq pc-2 in the central 20 pc (Ta- 
ble [6]), is among the highest compared to the mid- 
infrared survey of (ultra)lumin ous infrared galaxie s , star- 
burst s and Seyfert nuclei by iSoifer et al.r(|2000l I2OOIL 
120031 120041 ). In their 0'.'3-0'.'5 resolution survey, only 
three out of twenty-two nuclei were found to have a sur- 
face brightness (or its lower limit) at 10® Lq pc^^ or 
larger. All the three galaxies, Mrk 231, NGC 1275, and 
NGC 7469, have an AGN. Their mid-IR cores have (up- 
per limit) sizes of about 20 pc in two and 100 pc in 
one. For comparison, the luminosity surface densities 
of star formation-dominated nuclei are generally at or 
below 10"^ Lq pc"^. Between 10^ and 10® Lq pc"^ 
are ultraluminous infrared galaxies where both AGN 
and starburst m ay be hidden beh ind a large amount 
of gas and dust (ISoifer et al.l 1200 If) . More recent obser- 
vations by llmanishi et al.l ( 20 lit ) also support the mid- 
IR compactness and higher luminosity surface densities 
of AGNs compared to starbursts. There is a theoret- 
ical model explaining the apparent limit for starbursts 
around 10^ Lq pc~^. In it a nuclear gas disk opaque 
at far-IR will have that luminosity surface density when 
it supports itself by stellar radiation p ressure and self- 
regul ates to have Toomre's Q ~ 1 ([Thompson et al.l 
I2005D . This is for dust temper atures less than about 
200 K or for "warm starbursts" in lAndrews fc ThompsonI 
(|2011f) . For dust temperatures above 200 K, which the 
O'.'l core of NGC 4418 may have, and for a gas surface 
density of 10^ ® Mq pc~^ that corresponds to rggo /^m 1 
of the continuum core, the disk at the Eddington limit is 
in the "hot starburst" regime a nd will have a surface 
brigh tness of ~10®-^ Lq pc~^ ([Andrews fc ThompsonI 
I2OIII see their Equation 7). According to our data 
and this model, the nucleus of NGC 4418 has either an 
energetically-dominant AGN or the hot starburst. 

4.1.3. Constraint from Mass Budget 

Our estimates of the masses in the central 30 pc of 
NGC 4418 are summarized in Table [H The mass of 
young stars there, 10®-''='=°-^ Mq, is for a model young 
starburst having a luminosity of 10"'^^ Lq and a plausi- 
ble luminosity-to-mass ratio of lO^ o^O-S Mq"^ The 
mean star formation rate of this starburst would be 10- 
100 Mq yr~^ for a starburst age of 3 Myr. The mass of 
old stars in the central 30 pc, ~10^ '* Mq, is estimated for 
an average {^L*) Sa-type galaxy using the K-band pho- 
tometric parameters of .Graham fc Worlev. (,2008. their 



Table 3 ), their conversion form ulae to radial mass distri- 
bution ([Terzic fc Grahainll2005D , and the K-band mass to 
light ratio 1.06 Mq/Lqk for a 12 Gyr old stellar popula- 
tion with [Fe/II]=0 (Worthey 1994). From the viewpoint 
of mass budget, AGN dominance in luminosity would 
eliminate the need for young stars and help fit the sum 
of the component masses to the dynamical mass. The 
information currently at hand, however, does not let us 
conclude the AGN dominance from the mass budget be- 
cause the masses used here typically have an uncertainty 
of a factor of a few. The dynamical mass will be better 
constrained with higher resolution observations. 

4.1.4. Constraint from Time Scale 

The current tentative estimate of the outfiow kinemat- 
ical age, 10 Myr, is marginally longer than the age that 
an energetically-dominant starburst can have ( ij4.1.1l) . If 
this is firmly confirmed, it would be against the nascent 
starburst model. 

4.1.5. Constraint from Nh 

The gas shroud around the nucleus can shield the X- 
rays from the putative AGN because it is highly Comp- 
ton thick with Nn > 10^^ cm^^ to the galactic center 

(S121)- 

4.1.6. Provisional Verdict 

To summarize on the luminosity source, the large lu- 
minosity of the nucleus can certainly be due to a hidden 
AGN as far as the parameters L, L/M, I](Lboi), A'h, and 
mass budget of the nucleus are concerned. It would have 
an Eddington ratio ^0.3 with a factor of 2 uncertainty 
if the black hole and the bulge follow the normal scaling 
relation. This bolometric Eddington ratio is high com- 
pared to the median ratio of 10~^ among local type 1 
Seyferts (Ho 2008). In this sense the AGN mass accre- 
tion is rapid in NGC 4418 if the AGN is the dominant 
luminosity source. The black hole should be growing at 
a higher rate than in most low luminosity AGNs. The 
mass accretion rate to the black hole would be 0.1 Mq 
yj,-i jTQj. 20% radiative efficiency and the luminosity 
of IO^^Lq (i.e., Lboi ~ O.lMc^). The characteristic time 
scale for the black hole growth is 10® yr at this rate. This 
AGN gas consumption rate is much lower than the star 
formation rate required for the same luminosity because 
an AGN is more fuel-efficient. 

A compact young starburst that has a mass ^10® Mq, 
size < 20 pc, age - a few Myr, and L/M ~ 10^ Lq Mq"^ 
is an alternative that has not been excluded with our ob- 
servations. (The age could be up to 10 Myr as far as 
L/M is concerned, but starbursts older than 3 Myr are 
less favored because supernova explosions would be de- 
struct ive to the dusty core as argued by iRoussel et al.l 
()2003() .) While such a starburst would be at an edge of 
the parameter space allowed for starbursts, it is notewor- 
thy that the massive concentration of dense molecular 
gas is a favorable environment for star formation. This 
is true even if there is a major AGN. 

A combination of less luminous young starburst and a 
less luminous AGN is also possible as long as their total 
luminosity is about 10^^ Lq. It remains to be seen with 
more accurate mass, luminosity, and size measurements 
and the constraints mentioned above whether the nucleus 
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of NGC 4418 is in the parameter region where only a 
dominant AGN is ahowed. 

4.1.7. Consistency Check with Radio Data 

It is interesting to ask whether observations at radio 
wavelengths, which do not suffer from extinction by dust, 
are consistent with an energeticahy dominant starburst 
in the nucleus. As noted in §4.1.11 a starburst consis- 
tent with our SMA observations has a star formation 
rate of about 30 Mq yr~^ in the last 3.3 Myr so that 
the starburst population has a stellar mass of 10^ Mq, 
a luminosity to mass ratio of 10'^ Lq Mq~^, and the 
observed luminosity of Lboi ^ 10^^ Lq. A younger star- 
burst needs higher SFR to achie ve the luminosity. We 
note that the Lir-SFR re lation of lKennicuttl ()1998l ) and 
IKennicutt fc EvansI (poTl gives SFR w 15 Mq yr'^ for 
the W^^Lq. However, as one can see in Fig. [T7] (b), a 
continuous starburst at that SFR will be about 10 Myr 
old when achieving the luminosity and by that time the 
luminosity-to-mass ratio will drop by about 40% from its 
peak to ^600 Lq Mq~^ , assuming the standard Kroupa 
IMF. Thus, unless the IMF is truncated at low masses, 
our L/M constraint prefers higher SFRs than estimated 
from the standard Lir-SFR relation. 

We estimate a lower limit of SFR > 9 Mq yT~^ us- 
ing three constrai nts from the r a dio d ata and an SFR 
scaling relation of iMurphv et al.l ()2011[ ). The first con- 
straint is the high brightness temperature of the nu- 
cleus, >10^ K at ~2 GHz as reviewed in m.l[ at the 
compact nuclear source (~0'.'l-0'.'4). This exceeds the 
temperature of H II regions (~10'' K) and hence sug- 
gests the dominance of non-thermal synchrotron emis- 
sion. The second constraint is from the radio flux den- 
sities. The nucleus has 1.49 GHz (20 cm) and 4.86 GHz 
(6 cm) flux densities of 38.5 mJy and 26. 1 mJy, respec- 
tively , from a source w i th a 5. 0'.'4 size (jCondon et al.l 
Il99d iBaan fc Klockned I2006D : the resolutions of these 
observations are about 0'.'4 and 0'.'15, respectively. The 
20 cm luminosity of the nucleus is 5.3 x 10^^ W Hz~^ 
assuming isotropy of the radio emission. The third 
constraint is the spectral index between 20 and 6 cm, 
0^20-6 cm = —0.33 (for S oc i^"). It is closer to the spec- 
tral index of thermal free-free emission (a = —0.1) than 
to that of non-thermal emission (a « —0. 8) both from 
optically thin star-formin g regions (e.g., [Turner fc Hd 
Il994t IMurphv et al.l 1201 Ih . The apparent inconsistency 
between the first and third constraints can be reconciled 
if the radio emission is partly opaque, i.e., if either syn- 
chrotron self-absorption or free-free absorption by ther- 
mal plasma is significant. We therefore re-derived the 20 
cm flux density using the 6 cm data and a non-thermal 
spectral index of —0.8, while noting that this simplified 
correction will be insufficient if the 6 cm emission is also 
opaque. The corrected 20 cm flux density and luminos- 
ity are 67.2 mJy and 9.3 x 10^^ W Hz"-'^, respectively. 
The star formation rate of > 9 Mq yr~ ^ is from this 
20 cm luminosity and the formula (14) of IMurphv et al.l 
(|20Tll) . which in our case is (SFR/M0yr ^) = 6.64 x 
10-22(zy/GHz)-" **(L^/WHz-i). This formula was de- 
rived using the proportionality of supernova rate to the 
star formation rate (adopting the ratio between them 
from Starburst99) and the empirical relation between the 
supernova rate and the non-thermal radio luminosity. A 



similar star formation rate was derived in lC12l from radio 
data albeit without the opacity correction. 

The radio-derived SFR is smaller than needed for the 
starburst model at the face value. This may be for the 
following (at least) two reasons if the starburst model is 
correct. First, the radio SFR is a lower limit and can 
be higher as mentioned above. Second, the SFR may be 
underestimated because the starburst is very young and 
is still deficient of supernovae compared to older star- 
bursts for which the SFR formula above was derived. 
This nascent starburst is what Roussel et al. (200j|) ar- 
gued for from the relative weakness of radio emission 
compared to the infrared emission in this galaxy and 
by analyzing radio data in a way similar to ours except 
that they used lower-resolution data (beam > 45"). An 
obvious weakness of this second explanation is that the 
dominance of non-thermal emission (our first constraint) 
does not accord with the very young starburst. The ra- 
dio emission should be dominated by thermal emission 
from H II regions rather than by non-thermal emission 
from supernovae-induced cosmic rays as long as the nu- 
cleus is optically thin at radio wavelengths. One could 
still attribute the weak thermal emission to direct ab- 
sorption of ionizing photons by dust and to high-opacity 
of the free-free emission. Non-thermal synchrotron emis- 
sion can be less affected by these owing to the higher 
penetrating power of cosmic rays than ionizing photons. 
These qualitative explanations are not particularly at- 
tractive because of their complexity, but we do not yet 
have a proof to rule them out. 

In summary, the lower-limit SFR derived from radio 
data under the starburst-dominant model is about a fac- 
tor of three smaller than the lower-limit SFR needed 
for the luminosity of the nucleus. Although various 
uncertainties outlined above prevent us from excluding 
the starburst model, an energetically dominant starburst 
would need an elaborate model to explain the radio ob- 
servations. 

4.1.8. Implications for Previous Studies 

We also note implications of our new observations to 
a previous study of this energy source problem. Full 
reanalysis of all existing observations in light of our new 
d ata is beyond the scop e of this paper. 

Ilmanishi et al.l ()2004D used the presumed lack of shocks 
in the nucleus of this non-merging galaxy and the high 
HCN-to-HCO+ intensity ratio that they obtained toward 
the nucleus to infer that the warm H2 they observed in 
the nucleus is due to X-ray heating. They then esti- 
mated the AGN luminosity from the H2 line luminosity 
and reached a conclusion that the AGN is the dominant 
energy source in the nucleus. Our discovery of the red 
cone that we attribute to an out fiow from the nu c leus, as 
well as the discovery by iGonzalez- Alfonso et al.l ()2012al ) 
of the blue-shifted absorption presumably due to an in- 
flow, raise the possibility of shock contribution to the in- 
frared H2 lines. The HCN-to-HCO+ ratio does not give 
us the fraction of H2 luminosity due to an AGN. Regard- 
ing the HCN-to-HCO+ ratio itself, while its enhancement 
around AGN gained sup port from rnore re cent observa- 
tions at kpc resolutions (jKrips et al.ll2008l ). it has been 
also found to vary within the center of a single galaxy 
from one giant molecular cloud to another in the range 
covered by the low-resolution observations of starbursts 
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and AGNs (|Meier fc Turned I2012D . The latter authors 
suggested shock enhancement of HON citing stronger en- 
hancement of HCN tha n HCO^ in the molecular out flow 
from a protostar (Bachi iler fc Perez Gutierredll997[ ). In 
addition, for the particular case of NGC 4418, our detec- 
tion and analysis of rotational lin es fr om vibrationally 
excited HCN toward the nucleus (jSlOf ) further compli- 
cate the line ratio analysis by introducing the excita- 
tion path through the vibrationally excited state. At 
present, no estimate exists for the fractional contribu- 
tions of AGN X-rays and other mechanisms to the ob- 
served high HCN/HCO^ ratio. Thus these recent ob- 
servations make the assumed 100% contribution of X- 
ray heating to the H2 line luminosity uncertain. There- 
fore the AGN luminosity estimated from the H2 line 
luminosity is more u ncertain than it was at the time 
llmanishi et al.l ()2004l ) interpreted their data. This is not 
to say that the new observations are against the AGN hy- 
pothesis because the AGN dominance can probably still 
explain the observations and there are new observations 
that p refer an AGN for molecular ch emistry in the nu- 
cleus ([Gonzalez- Alfonso et al.ir2012bl ). As was the case 
in the radio analysis reviewed above, we have a problem 
that contributions of multiple (possible) mechanisms to 
observed properties make a conclusive data interpreta- 
tion diSicult. It seems that we still need to accumulate 
more constraints to obtain the breakdown of the lumi- 
nosity sources in the nucleus. 

4.2. The Gas/Dust Shroud of the Nucleus and 
Comparison with Other Observations 

The shroud of the 10^^ Lq source is a ~20 pc core 
having rgeo/jm ~ 1 and the surrounding >100 pc gas 
concentration that is generally warmer and denser to- 
ward the center. The Compton-thick core has a surface 
temperature ^150 K and a mean density '^5x10^ H2 
cm~^ f iJ3.3.2[) . The radial variation of gas properties is 
indicated by the emission size sequence of "CO > HCN, 
HCO^> dust continuum" fi )3.2p . the continuum proper- 
ties suggestive of "dense core + diffuse envelope" struc- 
ture (i j3.3.1l and l3.3.4l) . the larger velocity widths of lines 
with higher critical densities ( §3.4|) . and the dependence 
of line ratios and equivalent widths on area size f iJ3.5p . 
The same emission size se quence was observed in Arp 
220 ([Sakamoto et al.l [20091 ) . The declining temperature 
gradient from the center to larger radii has been an ex- 
planation for the deep 9.7 /im depression in the mid-IR 
spectrum ( §1.1|) . 

Our observations on the properties of the central gas 
concentration are generally consistent with those from 
the 1 mm observations in our companion paper (|C12f ). 
Leaving the comparison between the 1 mm and sub-mm 
observations to that paper, we only mention that th e de- 
convolved C0(2-l) peak brightness temperature in lC12l 
is 80 K to agree with the CO(3-2) peak Tb of 90 K and 
that '0121 also reports rad i al varia tion of gas properties. 

Gonzal ez- Alfonso et al.l ()2012aD made a detailed anal- 
ysis of the circumnuclear gas properties using Herschel 
spectroscopy at 50-200 /im. They detected dozens of 
high-excitation lines, mostly in absorption, from H2O, 
OH, O I, HCN, NH3, and their ^^O isotopologues at 100- 
300 km s~^ resolutions. Using these lines, they made a 
multi-component model in which the ISM is warmer in 
more compact components. Our observations are con- 



sistent with this trend. There are also agreements on 
the dust temperature around 150 K at r 10 pc and 
on the higher gas densities in more compact compo- 
nents down to the r 10 pc component (according to 
(■''dust, 200 /^m/size) in their Table 1). A difference between 
their model and our observations is that we detected 
about 0.2 Jy of 860 ^m continuum from a ~0'.'l (20 pc) 
core whereas their model predicts a quarter of that from 
a component of the same size and about the same tem- 
perature (150 K). Our larger flux density and brightness 
temperature lead us to a larger dust opacity for the core 
at 860 /im. In order for the 50-200 /im and 860 /im ob- 
servations to be consistent, it may be that the ISM in and 
around the core is not uniform to allow more leakage of 
emission from the warmer inner region than expected in 
the case of uniform shell s. Another possibility is t hat the 
dust /3 is smaller than 2. iGonzalez - Alfons o et al.l ()2012al ) 
also obtained an oxygen isotope abundance comparable 
to that in the solar neighborhood, [^'^O/^^O] ~500. As 
they noted, the lack of ^*0 enhancement through stellar 
processing is in accord with the luminosity source being 
an AGN and/or young pre-supernova starburst. 

4.3. Evolution of the Luminous Nucleus 

The formation and evolution of the the hidden lumi- 
nous nucleus must be closely linked to those of the mas- 
sive gas concentration at the nucleus. It is because the 
concentrated gas affects the luminosity source and its 
appearance through fueling and shielding. 

4.3.1. What drove the large amount of gas to the nucleus? 

The optical morphology of NGC 4418 (Fig. [1]) does 
not show significant disturbance or distortion to suggest 
a major merger. There is a blue companion VV655 at 
3' southeast (30 kpc on the sky) and at about the same 
redshift (2200 km s~^ j3- Although the two are probably 
physically associated with each other, the large separa- 
tion is unlike that in Arp 220 whose merger nuclei have a 
projected separation of only 0.3 kpc. For Arp 220 there is 
no doubt that the concentration of a large amount of gas 
to each nucleus was a consequence of the dynamical per- 
turbation of the major merger. In contrast, NGC 4418 
does not appear to have acquired its nuclear gas con- 
centration, and hence the nuclear luminosity, through a 
major merger although some properties of the nucleus 
are similar to those in Arp 220. 

A stellar bar and minor merger are milder mechanisms 
for gas concentration than major merger and galaxy in- 
teraction. For example, stellar bars are known to trans- 
port cold gas in galaxy disks, sometimes as much as >10* 
Mq, to the galactic centers (Sakamoto_el_al, 1999) . If 
the gas concentration in NGC 4418 were due to a stellar 
bar the situation may be similar to the one in the barred 
spiral galaxy NGC 3504 that has a single-peaked concen- 
tration of molecular gas with a scale length of 220 pc and 
a ma ss of 10^ Mq in the central 200 pc (Ke nnev et al.l 
Il993| ). The morphological classification of NGC 3504, 
(R)SAB(s)ab, is very close to that of NGC 4418 (Ta- 
ble [T]) including the same moderate bar-type distortion. 

^ What is cataloged in the atlas of interacting galaxies by 
IVorontsov-Velvaminovl 1)19771 ) is this companion itself and not the 
pair of VV655 and NGC 4418. 
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Currently our data of NGC 4418 show neither a kinemat- 
ical sign of bar-driven gas dynamics nor any clear sign of 
recent minor merger, although such signs may be found 
in the future in the somewhat disorderly gas motion in 
the central kpc (Fig. [7|). 

The nuclear gas concentration may also be because 
most of the gas came from bulge stars. The gas origi- 
nated from the galactic bulge through stellar winds and 
supernovae explosions would have little net angular mo- 
mentum. Therefore some of it will sink to the galactic 
center when cooled. A stellar population with a normal 
initial mass function returns about 30% of its mass to the 
interstellar space in 10^" yr pungwiert et al.ir200lD . The 
current bulge mass of NGC 4418 is about 0.6 x 10^° Mq, 
where we used the K band bulge luminosity in ij3.4l and 
the same luminosity to mass ratio as in ij4.1.3l The cumu- 
lative mass loss from the bulge stars should be therefore 
about 3 X 10^ Mq. Stars formed from this gas should 
have small specific angular momentum and hence fur- 
ther mass loss from them should keep the low specific 
angular momentum. The ~10^ AIq molecular gas cur- 
rently concentrated to the nucleus may be a fraction of 
the low angular mome ntum gas orig i nated t his way from 
the bu lge stars. fSee lYoung et al.l (|201lD . iDavis et"al] 
(|2011| ). and references therein for the general molecular 
gas content in local early type galaxies and discussion on 
their origins.) 

We note two more galaxies that have been recently 
found to host a compact and massive concentration of 
molecular gas at their centers. NGC 1377 is an iso- 
lated lenticular galaxy that shares infrared characteris- 
tics such as t he de ep 9.7 /xm absorption with NGC 4418. 
lAalto et all (poTl found a single-peaked '-^lOO pc-scale 
concentration of molecular gas at the nucleus and esti- 
mated a total gas mass in the central 400 pc of ~10^ Mq. 
NGC 1266 is another lenticular galaxy that has a strong 
concentration of molecular gas (~10*'^ -^0) in the cen- 
tral 100 pc without a clear sign of external dist urbance 
or bar-driven gas transport (|Alatalo et al.lf20Tll ). There 
may be a common mechanism to cause the compact mas- 
sive gas concentrations in these galaxies and NGC 4418. 
It may be the bulge-origin of the gas as mentioned above 
since these galaxies have large bulge fractions. 

4.3.2. Ongoing Evolution 

Our discovery of an outfiow feature from the center 
of NGC 4418 makes the nucleus to have both an out- 
flow and a previously known inflow. The inflow has 
been seen on our sightline near the galactic plane (i i3.4|) 
and the outflow is perpendicular to the galactic plane 
(at least at >500 pc scales). Such complex gas motion 
has been modeled around a galaxy nu cleus hosting both 
starburst and an a ctive nucleus ( e .g.. I\yada fc NormanI 
l200llWada|[20T^ . lOhsuga et aLl ()2005[ ) also simulated 
an outflow-|-inflow system at a smaller scale around a 
black hole and showed that the system can achieve the 
Eddington luminosity or higher. While we tentatively 
estimate the mass outflow ra te to be ^4 Mq yr~^, 
iGonzalez- Alfonso et al.l ()2012al) estimated an inflow rate 
of < I2M0 yr~^ from the redshifted absorption lines they 
found. The latter must be an overestimate because it as- 
sumed a spherically symmetric inflow without knowing 
the outflow. It is therefore possible that the net outflow 



rate is close to zero or even positive. 

The outflow from the deeply embedded nucleus may 
be driven at least at its root by the r adiation pres- 
sure on dust from the luminous nucleus (jScovi iig [2ool 
iMurrav et al.l[2005l ). The criterion for a radiation-driven 
outflow is given with the luminosity-to-mass ratio as 

Lhoi ^ 47rcG' 

Mdyn K 

where k is the mass opacity coefficient, c speed of light, 
and G the gravitational constant. For fully ionized hy- 
drogen plasma, k is the ratio of the Thomson cross 
section to the masses of proton and electron and is 
(Tt/("T'p + "^c) ~ 0.4 cm^ g^^- The minimum L/M 
for the outfiow is thus the Eddington limit, 3.3 x 10'* 
Lq Mq~ ^ . For 3 /im mass opacity coefficient of 50 
cm^ g~^ ()Scovillell2003[ ). where the opacity is due to dust, 
the minimum L/M for the outflow is 300 Lq Mq~'^. For 
50 /im where the SED of NGC 4418 peaks, the mass 
opacity coefficient is k = 2.5 cm g" -1 (IHildebrandlll98l 
and the minimum L/AI for the outfiow is 5000 Lq Mq~^. 
From these, a radiation-driven outfiow appears possible 
from the region where the peak of the SED is in near- 
IR; the peak shifts to longer wavelength further from the 
central source. Unlike the inner region where the radia- 
tion pressure may drive the outflow, the outer region of 
the outflow along the rotation axis of the galaxy may be 
driven by thermal pressure since the polar direction is 
the direction of the steepest pressure gradient. 

With the inflow/outflow rates of several Mq yi:~^ the 
ongoing evolution of the circumnuclear gas concentration 
is rapid in the sense that most of the ^10^ Mq molecular 
gas in the central 20 pc will be replaced in 10 Myr. Al- 
though the presumed accreting disk -I- outflow structure 
may stay longer with a continuous supply of gas from the 
larger gas reservoir outside the 20 pc core, the degree of 
nuclear obscuration may signiflcantly vary in the gas-flow 
time scale. The currently high L/Af can also change by 
an order of magnitude in 10 My for the starburst case 
and by a larger factor in much shorter time in the AGN 
case. These considerations suggest a rapid evolution and 
a short duration of the current state of the NGC 4418 
nucleus. 

Finally, NGC 4418 and the two galaxies mentioned in 
i )4.3.1l for their nuclear gas concentrations share an out- 
flow of molecular gas and/or dust from each nucleus. The 
nucl ear outflow from N GC 1377 is seen in CO kinemat- 
ics (|Aalto et al.ll2Q12D and, as in NGC 4418, as a red 
cone on a semi-minor a xi s in an optical c olor-i ndex image 
(|Heisler fc Vaded[l99l . lAlatalo et all (I20T1I) found the 
outflow of NGC 1266 from observations of CO, H I, and 
ionized gas and suggested it to be AGN-drivcn. There 
may be a common evolutionary relation between the 
compact nuclear gas concentrations and the outflows in 
these galaxies because a gas concentration is generally 
expected to feed nuclear activities that can drive an out- 
flow. 

5. SUMMARY AND CONCLUSIONS 

We made high-resolution 860 and 450 /im observa- 
tions of NGC 4418, a luminous infrared galaxy having 
a deeply embedded nucleus. Continuum and lines in- 
cluding C0(3-2), HCN(4-3), HC0+(4-3), and CS(7-6) 
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in the 350 GHz band were imaged at resolutions from 
O'.'S to 1"(=165 pc) in the central 34" (5.6 kpc). Contin- 
uum at 660 GHz (450 /im) was imaged at 0'.'2 resolution. 
We also examined SDSS optical images for circumnuclear 
features related to the luminous nucleus. Our key results 
are the following: 

1 . The interstellar medium traced by the submillimeter 
emission is found strongly concentrated toward the nu- 
cleus. CO (3-2) emission is detected only in the central 1 
kpc (6") of the galaxy with almost all the single-dish flux. 
More than half of the CO flux in the central kpc is in the 
central 165 pc (1"). Higher excitation lines other than 
CO are more concentrated toward the nucleus (FWHM 
~ 50 pc = 0'.'3), and the 860 fim continuum emission is 
even more compact. 

2. The peak brightness temperature of CO (3-2) in the 
nuclear gas concentration is as high as 90 K at 0'.'3 (50 pc) 
resolution, indicating the presence of hot molecular gas 
in the nucleus. Other indications of hot or highly excited 
molecular gas are lines from vibrationally excited HCN 
and HC3N with upper energy levels at 660-1070 K above 
the ground. 

3. The 860 /im continuum emission at the nucleus 
has a size of ^O'.'l (17 pc) in FWHM and deconvolved 
peak brightness temperature (Tb) of about 210 K for a 
Gaussian model. The diameter is 27 pc and Tb'^120 K 
for a uniform disk (sphere) model. The 860 /im opacity 
of this core is on the order of 1 and therefore the core 
is Compton thick (A'h ^ 10^^ cm~^). The bolometric 
luminosity of this dusty core is estimated from these pa- 
rameters to be iboi ^ 10^^" Lq. Most of the galaxy 
luminosity (10^^'^ Lq) is from this core. Its luminosity 
surface density is I](Lboi) ~ lO^-^iO-S Lq pc~^. 

4. Continuum emission at 450 /im also peaks at the 
nucleus but its flux density in our 30 pc («0'.'2) resolu- 
tion data is less than half of the total flux density from 
single-dish observations. A more extended region than 
the ~0'.'l core, presumably the region detected in C0(3- 
2), must have significant 450 /im emission. Extinction 
of the core by this colder envelope is a likely additional 
reason for the small core contribution at 450 /im. 

5. The dynamical mass within the central 30 pc is 
estimated to be ^2 x 10* Mq from the molecular gas 
velocities in the central 0'.'5. The nucleus thus has 
Lho\/ Mdyn ^ 500 Lq/Mq. Also, higher excitation lines 
are found broader than lower excited ones in the nucleus. 
This is likely due to faster gas motion and higher gas ex- 
citation parameters (e.g., density and temperature) to- 
ward the center. The CO line at r =0'.'25- 3" (40-500 
pc) shows complex velocity structure that likely contains 
radial motion in addition to rotation. 

6. The observed high L/M, S(i), and A'h as well as 
the nuclear mass budget are consistent with a hidden 
AGN to be the main luminosity source while they also 
allow a young starb urst. This confirms the proposition of 
iRoche et al.l (| 19861 ) and subsequent researchers. A dom- 
inant AGN would have a high Eddington ratio ^^0.3 for 
the black hole expected from the bulge luminosity. For 
star formation to achieve the observed L/M and S(Lboi), 
there must be a compact concentration of < 5 Myr old 
stars well within the central 20 pc with a total mass on 
the order of 10* Mq. Submillimeter data of higher qual- 
ity will further constrain the nature of the nucleus. 



7. We found in the optical data a U-shaped red feature 
along a semi-minor axis of the galaxy. We suggest it to 
be an outflow- induced cone seen in dust extinction (i.e., 
reddening) . We estimate the mass of the outflowing gas 
to be ~4 X 10^ Mq and the rate of mass outflow from the 
nucleus to be ~4 Mq yr~^, although the latter is tenta- 
tive. The luminous nucleus is blowing out its shroud in 
the polar direction while it is still fed from another direc- 
tion as suggested from the previously known redshifted 
absorption. Because these radial flows can replace most 
gas in the ^20 pc core in only 10^ yr, they should play a 
significant role in the evolution of the luminous nucleus. 
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TABLE 1 



NGC 4418 PARAMETERS 


Parameter 


Value 


note 


R.A. (J2000) 


12''26™54!612 


(1) 


Dec. (J2000) 


-00°52'39'.'41 


(1) 


Vsys [km s~^] 


2100 


(2) 


D [Mpc] 


34 


(3) 


Scale. 1" in pc 


165 




-t'8-lOOOMm [Lq] 


IQii.i 


(4) 




1010.3 


(5) 


Lb [Lqb] 


109-8 


(6) 


Hubble type 


(R')SAB(s)a 


(6) 


P.A. [°] 


60 


(7) 


i n 


62 


(8) 



Note. — (1) Position of the 860 
fim continuum nucleus measured in this 
work. The absolute positional uncer- 
tainty is estimated to be O'. 02 from the 
observations of a nearby quasar. Our 
pointing position was R.A. — 12^^26"^54^620 
Dcc.^-00°52'39''40 taken from the Sloan 
Digital Sky Survey Release 6 through the 
NASA Extragalactic Database. (2) Cen - 
troid velocity of CO(2-l) spectrum dSlOl) . 
Its formal uncertainty from the Gaussian 
fit is less than 1 km s""*" but see il3.4| 
(3) The galaxy distance we adopt. For 
V;ys(radio, LSR)^2100 km s"^, the galaxy 
re dshift corrected for l ocal flow in the way 
of iMould et al.| i2000D corresponds to an 
angular-size and luminosity distance of 33.5 
and 34. 1 Mpc, respectively, for Hq—73 

km S~^ MpC~^, Slmatter ^ 0.73, Slvacuum ^ 

0. 27. (4) From th e IRA S flux measurements 
in ISanders e t al. { 2003*) and the formula in 
Table 1 of [Slanders & Mir abeH 119961) . (5) 
Two Micron All Sky Survey, as listed in 
the NASA Extragal a ctic D atabase. (6) 
Ide Vaucouleurs et al.l <]199lh . (7) Position 
angle of the major axis from the Two Micron 
All Sky Surve y Extended Source Catalogue 
(•2MASS XSC: IJarrett et al.|2QO0l) . (8) Incli- 
nation of the galaxy disk. The minor to ma- 
jor axis ratio of b/a ^ 0.50 in 2MASS XSC 
is translated to inclination assuming the in- 
trinsic (edge-on) axial ratio of 0.2. 



TABLE 2 

Log of SMA observations 



No. 


UT date 


/lo 


A^ant 


array configuration 


-^baseline 


T225 


(Tsys) 


^obs 






[GHz] 




name 


pads 


[m] 




[K] 


[hr] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


1 


2009 Mar. 3 


348.742 


7 


VEX 


1,12,16,17,19,20,21 


40-509 


0.07 


237 


4.6 


2 


2010 Mar. 3 


347.698 


6 


VEX 


1,16,17,18,19,20 


85-509 


0.04 


153 


4.3 


2' 


2010 Mar. 3 


663.041 


5 


VEX 


1,17,18,19,20 


110-458 


0.04 


1228 


3.8 


3 


2010 Mar. 28 


347.698 


6 


sc 


1,2,3,4,5,6 


6-25 


0.05 


204 


2.1 


4 


2012 Feb. 5 


347.728 


7 


EXT 


1,9,12,14,15,16,17 


24-222 


0.04 


163 


4.9 



Note. — (3) Frequency of the local oscillator. Upper (lower) sideband is from 4.0 to 6.0 GHz above (below) 
this frequency for the first two tracks and from 4.0 to 8.0 GHz for the third and fourth tracks. (4) Number of 
available antennas. (5) SMA antenna conflguration. VEX— very extended, EXT— extended, SC— sub-compact. 
(6) Antenna locations. See lHo et al.| J2004i) for a map with the numeric keys. (7) Range of projected-length of 
baselines for NGC 4418. (8) Zenith opacity at 225 GHz measured at the Caltech Submillimeter Observatory 
adjacent to the SMA. (9) Median double sideband (DSB) system temperature toward NGC 4418. (10) Total 
integration time on the galaxy. 
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TABLE 3 
Data Properties 



ID 


Data Source 


SB. 


Emission 


Beam 


AV 


Noise rms. 










"x " 


km 


mjy beam ^ 


K(R-J) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


01 


1 


U 


HCO+(4-3) etc. 


0.34x0.24 


30 


27 


3.2 


02 


2 


u 


HCN(4-3) etc. 


0.36x0.26 


30 


20 


2.2 


03 


1+2 


L 


CO{3-2) 


0.35x0.26 


10 


22 


2.4 


04 


1+2 


D 


860 /im cont. 


0.35x0.26 


a 


2.4 


0.27 


05 


1+3+4 


U 


HCO+(4-3) 


0.67x0.45 


20 


15 


0.47 


06 


2+3+4 


U 


HCN(4-3) 


0.53x0.41 


20 


16 


0.72 


07 


1+2+3+4 


L 


CO{3-2) 


0.69x0.55 


10 


14 


0.38 


08 


3+4 


L 


CS{7-6) 


0.91x0.59 


20 


16 


0.32 


09 


3+4 


U 


4 GHz spectrum 


0.92x0.73 


30 


12 


0.17 


10 


3+4 


L 


4 GHz spectrum 


1.00x0.76 


30 


11 


0.15 


11 


2' 


D 


450 fim cont. 


0.23x0.15 


a 


65 


5.3 


Note. — Column 


(2): 


Data source. The numbers are those in Table [2] Column (3): 


Sideband. 



U— USB, L— LSB, D— double sideband— USB+LSB. Column (4): Main emission in the data. Column 
(5): Major and minor axes (FWHM) of the synthesized beam. Column (6): Velocity resolution. 
Columns (7) and (8): Noise r.m.s. intensity in the data. For lines this is for a channel of the width 
in (6). The noise in (8) is in kelvin and in Rayleigh- Jeans brightness temperature. 

^ The mean frequency and total bandwidth of the continuum are respectively 348.2 and 3.5 GHz for 
the 860 /im data and 663.0 and 3.9 GHz for the 450 fim data. 



TABLE 4 

Peak brightness temperatures and fluxes in the VEX data 



Emission 


max 


S, Sx (0'.'5 FWHM) 


/{0"5/6") 




K(R-J) 


K (Planck) 


Jy km s~^ , mJy 




(1) 


(2) 


(3) 


(4) 


(5) 


CO(J=3-2) 


82 


90 


240 + 3 


0.32 + 0.02 


HCN{J=4-3) 


24 


31 


80 + 3 


0.60 + 0.05 


HCO+(J=4-3) 


20 


28 


39 + 4 


0.48 ± 0.06 


Hi3cN(J=4-3) 


13 


20 


29 + 4 


0.71 + 0.13 


HC15n(J=4-3) 


8 


15 


17 + 4 


1.55 + 0.57 


HCN{«2=l^^,J=4-3) 


13 


20 


18 + 4 


0.60 + 0.15 


860 fim continuum 


17 


24 


168 + 3 


0.95 + 0.07 


450 fim continuum 


31 


45 


369 ± 37 




Note. — Columns (2) 


and (3): 


Peak brightness 


temperatures measured 


in the VEX 



data {data ID 01—04 and 11 in Table [3]l. No correction is made for missing flux and 
beam dilution. Temperatures calculated using the Rayleigh-Jeans approximation and 
the Planck function are given in (2) and (3), respectively. The maxTb is measured 
in the 10 km s~^ resolution data for CO and 30 km s~^ for other lines. Column (4): 
Line fluxes and continuum flux densities that we detected in the central O'/b (FWHM). 
They were measured from data ID 01—04 and 11 in Tablef3] after convolving them to the 
resolution. The line fluxes are integrated over 1900-2300 km s"^; for H^^CN and HC^^N 
we doubled the flux of half that width since the other half of the line is either outside 
of our bandwidth or blended with another line. Note that CO, HCN, and HCN('U2 = 1) 
lines are blended with a nearby HC3N transition. The la errors are for the thermal 
noise. We estimate additional flux-scaling error to be ~5 % for 345 GHz data and 40 % 
for 660 GHz data. Colu mn (5 ): Fraction of flux detected in our 0''5 data compared to 
the 6" resolution data in lSlOl . 
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TABLE 5 
Visibility Fitting Results 



source 


A, line 


flux density 


size 


reduced 




mm 


Jy 


arcsec 




(1) 


(2) 


(3) 


(4) 


(5) 


NGC 4418 


0.85 


0.205 ± 0.002 


0.101 ± 0.011 


0.63 


NGC 4418 


HCN{4-3) 


0.299 ± 0.006 


0.245 ± 0.012 


2.3 


NGC 4418 


HCO+(4-3) 


0.204 ± 0.006 


0.343 ± 0.028 


0.51 


Titan 


0.85 


3.41 ± 0.08 


0.848 ± 0.008 


1.9 


Titan 


0.88 


3.28 ± 0.06 


0.863 ± 0.007 


1.9 


J1222+042 


0.85 


0.55 ± 0.02 


0.049 ± 0.023 


2.8 


J1222+042 


0.88 


0.53 ± 0.01 


0.0 ± 0.0 


3.2 



Note. — Fitting of NGC 4418 used data from multiple array configurations; 
2010 VEX, SC, and 2012 EXT. The VEX data of HCO+ had too low S/N 
to be usable. Fitting of Titan and the quasar J1222+042 were to verify our 
data analysis. The size of Titan based on ephemcris is 0. 8425. Column (2): 
Continuum wavelength or line name. The 0.85 mm continuum is from USB 
and 0.88 mm is from LSB. Column (3): Flux density at zero baseline (i.e., 
total flux density). For lines, this is for the line emission averaged from 1900 to 
2300 km s~^. Column (4): The size is FWHM of a Gaussian for quasars and 
NGC 4418, and the diameter of a uniform-brightness disk for Titan. 



TABLE 6 

Luminosity Estimates of the Nucleus 



shape 


size 




Tb 


log Lboi 


logS(Lboi) 


logp(Lboi) 




mas 


pc 


K 




Lq pc-2 


Lq pc-3 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Gaussian 


101 ± 11 


17 ±2 


205 ± 47 


11.3 ±0.3 


9.0 ±0.4 


8.0 ±0.4 


Disk/Sphere 


162 ± 18 


27 ±3 


119 ±27 


10.8 ± 0.3 


8.1 ±0.4 


6.8 ± 0.4 



Note. — Column (1): Assumed shape of the source to calculate the peak brightness 
temperature. Disk is a circular disk (i.e., a sphere projected on the sky) with uniform 
surface brightness. Columns (2) and (3): Source size is the FWHM for the Gaussian 
model and the diameter for the Disk/Spliere model. Column (4): Brightness temperature 
calculated with the Planck function. This is the peak value for the Gaussian model and 
the value across the surface in the Disk/Sphere model. Column (5): Bolometric luminosity 
of the source calculated from the source size, temperature, and the Stefan-Boltzmann law. 
See ^3.3.31 Columns (6) and (7): The mean surface density and volume density of the 
luminosity within the size in the column (3). 



TABLE 7 
Line Velocities and Widths 



Line 



HCO+(4-3) 

HCN(4-3) 

CO{3-2) 

CS(7-6) 

N2H+(3-2) 

HCO+(3-2) 

HCN(3-2) 

CO(2-l) 



Vc 

km s 



2088.3 ± 4.2 
2085.8 ± 2.8 
2097.8 ± 0.2 

2088.7 ± 4.6 

2099.8 ± 5.4 

2088.9 ± 3.5 
2089.9 ± 3.0 

2100.4 ± 0.3 



FWHM 
km 



235.1 ± 9.8 

252.7 ± 6.6 

142.2 ± 0.5 

225.8 ± 10.8 

143.8 ± 12.8 

197.9 ± 8.2 
244.9 ± 7.9 
142.4 ± 0.7 



Note. — Centroid velocities and widths of 
molecular lines in the central 6" of NGC 4418. 
The velocities are in the radio definition and 
with respect to the LSR. These parameter s are 
from Gaussian fitting to the spectra in ISICH . 
Errors are ±lcr. HCN(4-3) is blended with a 
HC3N line on the low velocity side. 



SAKAMOTO et al. 



TABLE 8 

Mass Budget in the Central 30 pc 



type 


logM 


note 




[Mq] 




Mdy„(r < 15 pc) 


8.3 ±0.3 


(a) 




7.0 ±0.3 


(a) 


M^^l(r < 10 pc) 


~8.0 


(b) 


A^young stars (f < 10 pc) 


8.0 ±0.5 


(c) 


A^old stars (f < 15 pc) 


~7.4 


(d) 



Note. — Numbers with ~ may well have 
0.3—0.5 dcx uncertainties as others. 
^ Sec ^3A\ 

^ Sec 

^ See ^4.1.31 Most of the luminosity-bearing 
young stars must be well within the central 20 

See ^4. 1.31 For an average Sa galaxy with 
Mk ^ -23.8 mag. NGC 4418 is 0.3 times as 
luminous. 
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I 1 

3 kpc, 18" 



Fig. 1. — NGC 4418 optical image from the Sloan Digital Sky Survey (0.3-1 fim composite). The image size is 100" (16.5 kpc) on a 
side. North is up and east is to the left. 
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Fig. 2. — NGC 4418 spectra in the central 0'.'5 (80 pc; FWHM). Data from the VEX array configuration are used and continuum has 
been subtracted. Emission lines are marked at the redshift of 2100 km and six of them have velocity axes. The data are from two 
observations with off'set frequency coverages indicated with the two horizontal lines at the bottom. The data at the overlapped frequencies 
are averaged. The three error bars in each panel are for the three spectral sections (e.g., the middle one for the averaged section). 




0.0 ^ I cq o.q ^ 

340 341 342 343 352 353 354 355 

F(obs) [GHz] F(obs) [GHz] 



Fig. 3. — NGC 4418 spectra in the central ~1" (data ID 09 and 10 in Table[3|. No continuum subtraction has been made. The same 
spectrum is plotted twice for the LSB data to show weak features (left scale) and the full spectrum (right scale). Horizontal bars under 
strong lines indicate the velocity ranges of 1850-2350 km s-\ which are excluded when making continuum data. The central tick of each 
bar is at 2100 km s~^. 
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Fig. 4. — CO(3-2) channel maps made with the data from three array configurations. The synthesized beam of 0'.'69x0'.'55 (FWHM) is 
shown in the top-left panel. The contours start from ±35 mjy beam~^ = 0.95 K = 2.5 a and the nth contours are at 35ri^'^ mjy beam~^. 
Negative contours are dashed. The maximum intensity is 1.7 Jy beam~^ = 46 K at 2090 km s~^. The plus sign in each panel and the 
origin of the offset coordinates are at the position of the 860 fim continuum peak given in Table [T] 
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Fig. 5. — C0(3— 2) channel maps at the center of NGC 4418. Only the data from the VEX configuration are used. Contours are at —3, 
3, 6, 9, 12, . . ., 33(T, where 3ct = 65 mjy beam~^ = 7.3 K. Negative contours are dashed. The maximum intensity is 730 mjy beam~^ 
= 82 K at 2060 and 2070 km s'^. The plus sign in each panel and the origin of the offset coordinates are at the position of the 860 /im 
continuum peak given in Table[l] The synthesized beam of 0'.'35x0'.'26 (FWHM) is shown in the top-left panel. The last two panels show 
LSB continuum and the point spread function (i.e., CLEAN beam), respectively, in contour steps of 12.5% of each peak. 
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Fig. 6. — HCN(4-3) channel maps at the center of NGC 4418. Only the data from the VEX configuration are used. Contours are at 
(—1,1,2,3,4) X (50 mjy beam~^= 5.4 K = 2.5a). Negative contours are dashed. The maximum intensity is 23 K at 2020 and 2050 km s~^. 
The plus sign in each panel and the origin of the offset coordinates are at the position of the 860 /im continuum peak given in Table [Tl 
The synthesized beam of 0'.'36x0'.'26 (FWHM) is shown in the top-left panel. 
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Fig. 7. — CO(3-2) maps made using three array configurations, (a) Integrated intensity. The nth contour is at 60n^ K km s ^. The 
peak intensity is 9.4 X 10'' K km s~^. (b) Intensity-weighted mean velocity. Contour interval is 10 km s~^. The plus sign in each panel 
and the origin of the offset coordinates are at the position of the 860 fira continuum peak given in Table [l] The black ellipses show the 
FWHM of the synthesized beam. 
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Fig. 8. — Submillimeter line and continuum images at the center of NGC 4418; (a)-(d) are CO(3-2), (e)-(h) are HCN(4-3), (i) is 
HCO+(4-3), (j) is HCN(t)2=l^-^,J=4-3), and (k) and (1) are 860 and 450 fiui (DSB) continuum, respectively. Data from the VEX 
configuration are used. The position of the 860 fim continuum peak, which is given in Table [T] is marked with the plus sign in each panel 
and is used as the origin of the offset coordinates. The black ellipse in each panel shows the FWHM of the synthesized beam. The images 
labeled with 'momO', 'moml', 'mom2' show integrated intensity (i.e., 0-th moment) in K km s~^, intensity-weighted mean velocity (1st 
moment) in km s~^, and velocity dispersion (2nd moment) in km s~^, respectively, and those labeled with 'max T^,' show peak brightness 
temperature in K. The contour information is as follows, (a) The nth contour is at 0.5n^'* X 10^ K km s^^. The peak value is 18 X 10^ K 
km s~^. (b) 10 km s~^ steps. The contour crossing the continuum peak is at 2105 km . (c) 10 km s^^ steps, (d) 10, 20, 40, 60, and 80 
K (lo-=2.4 K). The peak is 82 K. (e) 0.78n 10^ K km s"!. The peak is 6.7 X 10^ K km s"!. (f) 25 km s"! steps. The contour crossing the 
continuum peak is 2075 km s"!. (g) 10 km s'^ steps, (h) (1, 2, 3) X (6.6 K = 3 ct). The peak is 23 K. (i) (1, 2, 3) X 0.76 X 10^ K km s'^. 
(j) The contour is at 0.74 X lO"' K km s^^. (k) 0.8n^-^ K, starting at Scr. The intensity peak is 17 K. (1) Contours are at (—3,3,5) X (Ict 
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Fig. 9. — Visibility fitting results, (a)-(c): Continuum data of NGC 4418, J1222-I-042 (a test point source), and Titan (a source with 
known size). The latter two show 2010 VEX data alone. Dotted lines are Gaussian fits for (a) and (b) and a fit with a circular disk of 
known diameter for (c). (d)-(e): Line data of NGC 4418 integrated over 1900-2300 km s"^ after continuum subtraction. Data shown in 
blue, green, and red are from the SC, EXT, and VEX configurations of the SMA, respectively. Dotted lines are models fitted to the data. 
See Table [5] for the fit parameters. The CO data could not be reasonably fit with a single Gaussian. 
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Fig. 10. — Spectral energy distribution of N GC 4418 in its Rayle igh- Jeans part. Filled data points are from SMA observations; blue ones 
are from the short baseline observations in .Sakamoto et aP lj201CI ) and red ones are from long baseline observations. Error bars are ztlcr. 
The dashed line is a power law fit to the SMA data between 200 and 400 GHz. It has a spectral i ndex of a = 2 .55 it 0.18 (for 5^ oc u") and 
is extrapolated to a wider frequency range to guide eyes. Non-SMA data are from the literature fBrauher et al. 2008; Sanders et al. 200^; 
I Yang fc Phillipsll2007l ; IHoche fc Chandlen[T993 : IDunne et al.ll200CII ; llmanishi et al.1120041 '1. No subtraction for line contamination is made for 
the bolometer data between 100 GHz and 1 THz. 
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Fig. 11. — (Top) Sampling of our 660 GHz observations in the u— u plane. (Bottom) Fraction of &ax recovered in the simulations of our 
observations for Gaussian sources with various sizes. 
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Fig. 12. — Integrated intensity and mean velocity maps of HCN(4-3), HCO+(4-3), and CS(7-6) emission. Data from tliree array 
configurations are used except for CS that was observed in two. The nth contour is at 70n^'^, 50n^''', and 30ra^'^ K Icm in (a), (c), and 
(e), respectively. Velocity contours are in 10 km steps. The position of the 860 fim continuum peak, whose coordinates are in Table [Tl 
is marked with the plus sign in each panel and is used as the origin of the offset coordinates. The black ellipses show the FWHM sizes of 
the synthesized beams. 
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Fig. 13. — Position-velocity diagrams of molecular lines at the center of NGC 4418. Each PV cut is through the continuum peak and 
along the position angle of 45°. The origin of the abscissa is the continuum peak. The black rectangles at the bottom-right show the 
spatial and velocity resolutions of the data. Contour steps are 8 K (3.3(t) for CO, 3.5 K (1.6(t) for HCN, and 5 K (1.6(t) for HCO"*". The 
highest- velocity emission in the CO diagram is from H^^CN, and the low-velocity signal in the HCN diagram likely contains HC3N(39-38). 
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Fig. 14. — (top) Line centroid velocity (Vc) and FWHM for the lines in TableO (bottom) Critical density (ricrit) at 50 K and FWHM 
for the same lines. Error bars are ±4(7 for CO and ±lcr for others. Molecular parameters for the ricrit are from the Leiden database 
HSchoier et aU\^OB) . 
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Fig. 15. — CO(3-2) and HCN(4— 3) line profiles in the central arcsecond of NGC 4418, sampled on a grid of 0'.'12 spacing. The origin 
of the offset coordinates is at the 860 fim continuum peak, the spectrum observed at which is labeled with C. Velocity resolution of the 
spectra is 10 and 30 km s"'^ for CO and HCN, respectively. The rms noise per channel is 2.4 and 2.2 K for CO and HCN, respectively. 
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Fig. 16. — Distribution of optical color in the center of NGC 4418. Flux ratio involving four SDSS bands, (i' + z')/{g' + r'), is plotted 
in such a way that regions brighter in the longer wavelengths i' and z' than the shorter wavelengths g' and r' are shown in red while low 
ratio regions are shown in blue. The resolution of this image is I'.'S. The origin of the offset coordinates is at the 860 /xm continuum peak 
and is marked with a plus sign. North is up and east is to the left. 
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Fig. 17. — Evolution of starburst parameters simulated using Starburst99 and the IMF of lKroupal 1200^ '). which has a power law index 
of —1.3 for 0.1-0.5 A/0 and —2.3 for 0.5-100 Mq. (a) Luminosity-to-mass ratios for single-age (i.e., instantaneous) starburst populations 
with three different IMF mass ranges. Also plotted is the Eddington L/M for fully ionized hydrogen plasma; the mass is the black hole 
mass in the context of AGN. (b) Luminosity and luminosity- to-mass ratio of continuous starbursts with different SFRs. Thin solid lines 
are starburst luminosities (left axis) for SFRs between 1 and 300 Mq yr~^. The thin dashed line is the starburst luminosity for SFR=10 
Mq yr~^ and for the Kroupa IMF truncated to the range of 1—100 Mq. The thick lines are the luminosity-to-mass ratios (right axis); the 
dashed one is for the same truncated IMF as above. The L/M ratio is independent of the SFR. 



